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Abstract 

The curvelet transform is a directional wavelet transform over M", which is used to analyze 
functions that have singularities along smooth surfaces (Candes and Donoho, 2002). I demon- 
strate how this can lead to new quantum algorithms. I give an efficient implementation of a 
quantum curvelet transform, together with two applications: a single-shot measurement proce- 
dure for approximately finding the center of a ball in K", given a quantum-sample over the ball; 
and, a quantum algorithm for finding the center of a radial function over M", given oracle access 
to the function. I conjecture that these algorithms succeed with constant probability, using one 
quantum-sample and 0(1) oracle queries, respectively, independent of the dimension n — this 
can be interpreted as a quantum speed-up. To support this conjecture, I prove rigorous bounds 
on the distribution of probability mass for the continuous curvelet transform. This shows that 
the above algorithms work in an idealized "continuous" model. 

1 Introduction 

One of the most remarkable demonstrations of the power of a quantum computer is Shor's algorithm 
for factoring and discrete logarithms [26] . This has motivated many researchers to try to generalize its 
key components — the quantum Fourier transform over Z^v, and the algorithm for period-finding — to 
solve other problems [13j. In particular, this motivated the study of the quantum Fourier transform 
and the hidden subgroup problem (HSP) on non-Abelian groups, as a route to solving certain lattice 
problems and the graph isomorphism problem [23^ El [25] . 

In this paper we study a different generalization of the Fourier transform, namely the curvelet 
transform over M" [lOj. This is a kind of "directional" wavelet transform, which can resolve features 
over the spatial and frequency domains simultaneously. A curvelet basis function resembles a wave- 
packet, with high-frequency oscillations in one direction (like a plane wave e^^'^, as in the Fourier 
transform on M"), but which is also supported on a small region of space (unlike the plane wave). 
We show that this leads to fast quantum algorithms for some new classes of problems, outside the 
framework of the HSP. To the best of our knowledge, this is the first attempt to design quantum 
algorithms based on the curvelet transform. 

Intuitively, the curvelet transform is helpful in analyzing functions on M" that are discontinuous 
along (n — l)-dimensional surfaces. If a function / is discontinuous along a surface S, then its curvelet 
transform Tf will be "large" at those locations (b,9), where 6 is a point on S and 6 is the vector 
normal to S at b. The set of all such pairs (6, 6) is called the "wavefront set" of /. 
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The basic model for a quantum algorithm using the curvelet transform is as follows: first prepare 
a quantum state Ylx^M" which is a weighted superposition of points in M"; then apply the 

quantum curvelet transform, to get the state X^gg rj(6, 9)\b, 9); and finally measure b and 9. For the 
time being, we ignore implementation issues, such as how to discretize and how to compute the 
quantum curvelet transform efficiently. The more basic question is whether we can find functions / 
such that we can prepare the initial state efficiently, and such that a measurement in the "curvelet 
basis" yields useful information. 

One example consists of letting / be the indicator function of a ball in R". We can efficiently 
prepare a uniform superposition over points in a ball, using the techniques of [2 [17]. (This is called 
"quantum-sampling," and it can be done more generally, e.g., for convex bodies.) Then, measuring 
in the curvelet basis extracts information about the center of the ball. 

Another example consists of choosing / to be the indicator function of a spherical shell in M". 
This is motivated by the problem of finding the center of a radial function on M". Let G be a radial 
function, centered around some unknown point c. Prepare a uniform superposition over a large 
region of space, then compute G{x) and measure it; this produces a uniform superposition over one 
of the level sets of G, which is a spherical shell centered at c. Then, measuring in the curvelet basis 
extracts information about the location of c. 

The goal of this paper is to make this intuition precise. We can interpret Fj as a wavefunction, 
and we want to show that its probability mass |Fjp is concentrated near the wavefront set. We 
can prove this for the continuous curvelet transform, for the two cases of interest, where / is the 
indicator function of a ball or a spherical shell in M". In these cases, (6, 9) is near the wavefront set 
with high probability. This implies that the line {b + X9 \ X € M} passes near the center of the ball 
or spherical shell. 

Next, we give an efficient implementation of a quantum curvelet transform. (This is a dis- 
crete version of the transform described above, acting on superposition states.) Then we propose 
polynomial-time quantum algorithms for two problems: (1) given a single quantum-sample over a 
ball in M", find the center of the ball, with accuracy it A where A is a constant fraction of the radius 
of the ball; (2) given oracle access to a radial function / that is centered around some unknown 
point c G M", find the point c exactly (i.e., with accuracy itA in time poly(log-^), assuming that 
the function / fluctuates on sufficiently small scales). 

For the first problem, we conjecture that our quantum procedure succeeds with constant prob- 
ability, while the best classical procedure succeeds with probability that is exponentially small in 
n. Classically, this problem is hard because the volume of a ball in R" is concentrated near its 
surface. But this same fact is helpful to the quantum curvelet transform, which works by finding a 
line normal to the surface of the ball. 

For the second problem, we conjecture that our quantum algorithm uses only a constant number 
of queries, whereas any classical algorithm requires J7(n) queries. Intuitively, this is because the 
curvelet transform uses constructive interference to find a direction in M" from just one query. 

We then prove that these algorithms work in an idealized "continuous" model — this follows 
from our rigorous results on the continuous curvelet transform. However, we do not have a rigorous 
analysis of the effects caused by discretization, though we can argue that these should be small. 

These examples demonstrate that one can use the curvelet transform to obtain a quantum speed- 
up. These examples are artificially simple, in order to allow a rigorous analysis. But the underlying 
idea — using the curvelet transform to find normal vectors to a surface — should work on more com- 
plicated geometric objects. 

1.1 Technical Contributions: First, in Section 2, we define the continuous curvelet transform over 
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M". This generahzes the definition over given in [8j. Given a function f{x), the continuous curvelet 
transform returns a function T f{a,b,9). Here, x € M" represents a "location," while < a < 1 is a 
"scale" (smaller values denote finer scales, larger values denote coarser scales), b E M" is a "location," 
and 9 S S^~^ (the unit sphere in R") is a "direction." 

Next, we study the distribution of probability mass \Tf\'^ over different values of (a, 6, 9). This is 
technically quite difficult. ^f{a, b, 9) is defined by an oscillatory integral, and while there are various 
methods for bounding the asymptotic decay rates of such quantities [28l [7] , we need non-asymptotic 
bounds on the total probability mass in a given region. In Section 3 we develop some tools for proving 
such bounds, in the case where / is a radial function. Then, in Sections 4 and 5, we specialize to 
the case where / is the indicator function of a ball or a spherical shell. Here, the analysis relies on 
powerful classical results that bound the oscillation and decay of Bessel functions [T| I30j. 

In Section 4 we let / be the indicator function of a ball in R", with radius centered at the 
origin. We expect that, after applying the curvelet transform, b and 9 will be concentrated near the 
wavefront set of /: that is, b will be concentrated near the line {X9 \ A G R}, at about distance /3 from 
the origin. Furthermore, we expect that b will become more tightly concentrated, the smaller the 
value of a. We show that this is essentially what happens. In particular, with constant probability, 
the distance from b to the line {A^ | A € R} will be at most a constant fraction of /3; remarkably, 
this holds independent of the dimension n. 

In Section 5, we let / be supported on a thin spherical shell, having radius /? and thickness 
(5 <C /3. Here, after applying the curvelet transform, we get a qualitatively similar behavior of a, b 
and 9. Quantitatively, however, we find that we can observe much smaller scales a, on the order of 
5/P; and thus we can find the center of the shell with much greater precision. Essentially, by making 
the shell extremely thin, we can find its center with arbitrarily high precision. 

Finally, we turn to the discrete curvelet transform, and quantum algorithms. In Section 6 we 
give an efficient implementation of a quantum curvelet transform. This uses ideas from the fast 
classical curvelet transform [6j. However, there is a new complication in the quantum case: we need 
to prepare certain states which are superpositions of different scales and directions (a, 9). We design 
families of curvelets that allow this step to be performed efficiently, and that have similar analytic 
properties to the curvelets used in Sections 3-5. 

In Section 7, we formally define the two problems mentioned earlier: estimating the center of a 
ball, given a single quantum-sample state; and finding the center of a radial function, given oracle 
access. We present quantum algorithms for these problems, and use our results from Sections 4 and 
5 to prove that the algorithms work in a continuous model. We also sketch a classical lower bound 
for finding the center of a radial function. 

This paper omits most of the proofs, due to lack of space. The proofs can be found in the full 
version [21]. Also, note that this paper contains some additional results there were not present in 
the first version of [21] . This paper contains an improved algorithm for finding the center of a radial 
function, and a classical lower bound for that problem. 

1.2 Related Work: Curvelets over R^ and R^ have been studied as a tool for image processing 
and simulating wave propagation [10^ [9l [HI |20] . The curvelet transform is also related to older ideas 
from harmonic analysis; see, e.g.. Smith [27]. It can be viewed as an algorithmic implementation of 
a technique known as second dyadic decomposition [28j . 

There are a few rigorous results on the behavior of the curvelet transform which are similar in 
spirit to our work [lOl El |5] . These results apply to much broader classes of functions, but they are 
only known to hold over R2 (or R3 in some cases) . Although one would expect them to generalize in 
some fashion to R", it is perhaps surprising that the scaling with n is as favorable as we find here. 
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In connection with quantum algorithms, there has been some work on the quantum wavelet 
transform [3H [T5\ [321 [T6] . But the curvelet transform on R" is quite different from the "ordinary" 
wavelet transform on M^, which consists of a product of 1-D transforms. The ordinary wavelet 
transform on can detect the locations of discontinuities, but it cannot resolve directions. 

The geometric problems studied in this paper are reminiscent of some recent work on finding 
hidden nonlinear structures, although the details are different. Shifted subset problems [12\ [22] 
involve translational invariance, so the natural tool for solving them is the Fourier transform, rather 
than curvelets. Hidden polynomial problems jl2l [H] resemble our problem of finding the center of 
a radial function. However, they are much more general (and thus much harder), and they are set 
over a finite field rather than W\ 

These problems can also be studied from the perspective of quantum state discrimination [llj . 
e.g., what is the optimal quantum measurement for estimating the center of a ball? However, we 
stress that our algorithms using the quantum curvelet transform are computationally efficient. 

Finally, we recently became aware of a quantum algorithm for estimating the gradient of a 
function / on M", using only 0(1) queries [19] . This is quite similar to our algorithm for finding the 
center of a radial function — it is like applying the curvelet transform at a single location b. Viewing 
this as a curvelet transform has the advantage of providing a more general framework, where one 
can do this procedure on an arbitrary input state. Also, note that our emphasis in this paper is on 
functions / that are not smooth — in this case, the gradient is not well-defined. 

2 The Continuous Curvelet Transform 

We begin by defining the continuous curvelet transform over M"-. This generalizes the definition of 
[8j over M^. Given a function /(x), the continuous curvelet transform returns a function T f{a,b,6). 
Here, x £ represents a "location," while < a < 1 is a "scale" (smaller values denote finer scales, 
larger values denote coarser scales), b £ M" is a "location," and £ S^~^ (the unit sphere in W^) is 
a "direction." All functions return values in C. 

Intuitively, the curvelet transform decomposes / into pieces corresponding to different scales a 
and directions 6; we can view ^f{a, b, 9) as a family of functions, indexed by a and 6, each representing 
some "piece" of f{x) (note that the variables b and x both represent locations in space). To get one 
such piece of /, we will take its Fourier transform, and then multiply by a "window function" q 
which is defined over the frequency domain. 

More precisely, the curvelet transform is defined to be 

Tf{a,b,9):= [ f{k)x,^<k)e^^''-'dk. (1) 

Here, / is the Fourier transform of /, and g is a function that is smooth, real, non-negative, and 
supported on a "sector" of frequency space S^^C M". We will describe these sectors below. Before 
doing so, we remark that the curvelet transform consists of (1) taking the Fourier transform of /, (2) 
separating / into pieces corresponding to different scales and directions, and (3) taking the inverse 
Fourier transform. This description suggests how to compute the curvelet transform efficiently. 
The "sector" ^ is roughly given by the intersection of the cone centered around the vector 

9 with angular width ^/a, and the annulus with inner radius 1/a and outer radius 2/a. Thus, the 
"piece" of / at scale a and direction 9 is somewhat like the restriction of / to frequencies k ~ {\/a)9 
(which represent oscillations in direction 0, at higher frequencies when the scale a is small). Note, 
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however, that the sector has dmicnsions 1/a x , so its shape is not constant — the sector 

becomes longer and narrower when the scale a is small. 

This construction can also be understood from a second perspective. We can define a family of 
curvelet basis functions l^-^§ StS follows. The curvelet at location h is obtained by translation from 

the curvelet at location 0, that is, 7^5^(2?") := 7^ e^'-^ ~ '^^^ curvclct at location is defined in 

terms of its Fourier transform, which is simply the window function ^, that is, l^^gik) := ^(^)- 
It is easy to check that the curvelet transform defined earlier is equivalent to taking inner products 
with this family of curvelet basis functions: Vf{a,h,6) = ilabe^^^' 

Now we can see how our choice of the window function ^ implies (and is motivated by) certain 
properties of the curvelet basis functions 7^ ^. Since x„ g is smooth, the "f^j^ff are rapidly decaying. 

Also, each 7^ has high-frequency oscillations in the 9 direction, and is essentially supported on a 

plate-like region, centered at location b, orthogonal to 9, with dimensions a x {^/a)"'~^. Intuitively, 
7^ g ^ resembles a plane-wave in direction 9, localized around the point b. 

Finally, we define the window function ^ as follows. Write A; using spherical coordinates 

{r,(j)i, . . . centered around the direction 9, so that (pi is the angle between k and 9. Then 

define Xaff^^) '■— ^(^or)l^(0i/v^)Aa(</'i). Here A is a constant, which can be chosen freely; we 
will explain how to set it later. W is a radial window function, real, nonnegativc, supported on 
the interval [1/e, 1], and satisfying the admissibility condition J^W{r)'^^ = 1. F is an angular 
window function, real, nonnegativc, supported on the interval [0, tt/2], and satisfying the admissibility 
condition Jg„-i V{(pi)^da{(pi, . . . , (pn~i) = ^, where da denotes integration over the unit sphere S'^~^ 
in M". ^ ^^^^ 

Afl is a normalization and adjustment factor: Aa{(pi) := q("-+i)/^ . This is 
needed because the volume of the sector ^ (on which x^ is supported) changes with a. Note 

that, when (pi is small, Aa(^i) ~ a^""''^^/'^. This is the main point where defining curvelets over is 
more complicated than over M?; note that in dimension n = 2, Aa(<pi) = 0*^"+^^/^ = a^/^ exactly. We 
remark that a simpler approach would be to use a constant normalization factor that only depends 
on a and not (pi; however, our more complicated construction will be more convenient in the later 
sections of this paper. 

We now state some basic properties of the curvelet transform. First, note that the curvelet 
transform works primarily on the high-frequency components of /, which correspond to fine-scale 
elements (a small). The constant factor A, mentioned above, sets the low-frequency cutoff value, 
which corresponds to the coarsest scale (a = 1). For convenience, here we assume that / has 
no low-frequency components below the cutoff value. In practice, when / has such low-frequency 
components, the curvelet transform leaves them unchanged, and simply returns them as a residual 
function /res- 

Next, we define the reference measure djj,{a,b,9) := ^t+T^^^ ■ This weights the contributions 
of Tf{a,b,9) differently according to the scale a. Intuitively, this is needed because the sectors 
S^^ do not cover the frequency domain uniformly, and the translations of a curvelet ^^^^g (for 

fixed a and 9) to different locations b do not cover the spatial domain uniformly. Rather, (a, 6, 9) 
should be "sampled" in a certain way. To see this, write dfj,{a,b,9) = ^ ^(n+^i) /2 atn-i)/2 ■ Note that 
^ = d(loga), suggesting that we should sample log(a) at uniform intervals, i.e., we should set a 
equal to powers of 2; we should sample 6 on a grid in whose cells have size a x {■\/a)'^~^\ and 
we should sample ^ on a mesh on S""^ whose cells have size (y^)"~^. Later, when we construct 
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the discrete curvelet transform, we will use this sampling trick for a and 9, in place of the reference 
measure. 

Then we have the following theorems: 

Theorem 1 Suppose that f{k) = for all \k\ < 1/A. Then we can recover f from its curvelet 
transform T/; /(f) = /^^-^ Tf{a, b, ^)7^ gg-(f )(i/i(a, b, 6). 

Theorem 2 Suppose that f{k) = for all \k\ < 1/A. Then the curvelet transform preserves the 
norm: \f{x)\'^dx = J^^^ \r f{a,b,9)\'^dn{a,b,e). 

These are straightforward generalizations (to the case of M") of results in [8j. We sketch the proofs 
in Appendix A. 

3 The Curvelet Transform of a Radial Function 

\T f{a,b,6)\'^dn{a,b,9) can be interpreted as a probability density over the different scales, locations 
and directions {a,b,6). In this section we will develop some tools for understanding where this 
probability mass is concentrated. We will consider the case where / has rotational symmetry. Though 
there is no simple analytic expression for Fj, we can deduce certain properties from symmetry, and 
we can upper-bound the variance of b (this latter point is our main result). 

Let / be a radial function, f{x) = /o(|x|). Its Fourier transform is also radial, f{k) = Fo(|/c|), 
where Fo{p) = ^(„^^)/2 J{n-2)/2{'^'^P'^)ff){^)'''^^'^dr^ and J is a Bessel function (see, e.g., l29j). We 

assume that / is normalized so that \f{x)\^dx = 1. 

When / is radial, Tf has the following symmetries: Tf{a,b,6) = T f{a, —b, —6), and for any 

rotation R, Tf{a,b,9) = T f{a, R{b), R{e)). 

We make a particular choice for the radial and angular windows W and V. These windows are 
smooth, which is necessary in our analysis of the variance of b. We let W{r) = [Cw sin(7rlogr)^ 

if 1/e < r < 1; otherwise] and V{t) = [Cy cos{t)'^ if < t < 7r/2; otherwise], where = ^/8/3, 

and a = ^/^^^■ 

3.1 The probability of observing a scale a: First, we claim that the probability of observing 
a fine-scale element a < rj is essentially given by the amount of probability mass of / at frequencies 
above l/(Ar/): Pr[a < 77] = \T f{a,b,e)\'^dfi{a,b,e) > fik\>i/{Xn) \fik)\'^dk. This follows from the 
same argument used to prove Theorem [2l In the case of a radial function, we write this as: 

POO 

Pr[a < ??] > 5o / Foipfp^-'dp, (2) 

Jl/{Xri) 

where 5*0 = ^{71/2) surface area of the sphere S"'~^ C M". 

3.2 The location b and direction 9: We claim that the location b has expectation value 0. To 
see this, observe that E(6j) = — E(6j), due to the reflection symmetry of Fj; thus E(6j) = 0. Note 
that this remains true when we condition on the value of a. 

Also, we claim that the direction 9 is uniformly distributed. This follows from the rotational 
symmetry of Fj. Note that this remains true when we condition on the value of a, and when we 
condition on the value of b ■ 9 (since this preserves the rotational symmetry) . 
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3.3 The variance of 6 perpendicular to 6: Finally, we seek to upper-bound the variance of b, in 
the directions perpendicular to 9, as well as parallel to 9. These results are rather complicated, so 
we defer most of the details to Appendix B.l. However, these results are a basic component of our 
proofs in Sections 4 and 5, so we will sketch some of the calculations. 
First, the variance of b perpendicular to 9 is: 

Eib^{I-99^)b)= [ C f {b^{I-99^)b)\Tf{a,bMdb^da{9). (3) 

Note that a similar formula holds when we condition on observing a < rj. 

We can take advantage of rotational symmetry to do the 9 integral. Fix a vector u = (1, 0, . . . , 0), 
and for each 9, let i? be a rotation that maps 9 to u. Then we can replace the expression inside 
the integral with {R{b)^{I — uu^)R{b))\T f{a, R{b),u)\'^. Then change variables b ^ R~^(b). The 
integrand is now independent of 9, so we can do the 9 integral. We get: 



E{b^{I-99^)b) = So l' I {b^{I-uu'^)b)\Tf{a,b,u)\^db^. (4) 

JO JK" " 

Now the key idea is to replace integration over the spatial domain with integration over the 
frequency domain, via Plancherel's theorem. (Recall that curvelets are defined more simply over the 
frequency domain.) We introduce some new notation, ^^^(b) := Tf(a,b,9). By equation ([1]), the 

Fourier transform of g is given by g{k) = f{k)x^ gi^)- 

Let Ik denote the innermost integral in equation ([4]). Then Ik = Y2^=2 Ir^ \bj\'^\^a,u{b)\'^db. 
Using Plancherel's theorem, and symmetry with respect to rotations around the u axis, we can write 

We can expand out the integral on the right hand side, as follows. Using spherical coordinates 
k = (r, (/>!,... , (pn-i), we write ^a,u{k) as a product of a radial part and an angular part: ^a,u{k) = 
L{r)M{(j)i), where L{r) = Fo{r)W{\ar), and M{(j)i) = V{(j)l/^/a)Aa{c|}l). Then we have 

i-^^adk) = L'{r)M{^i)^^ + L{r)M\<p^)^, (5) 

where ^ = sin^i cos(/>2, and = S2^Jh^2lil^ go we get: (note that ^a,u{k) is real) 

= TTTv? / [L'{r)M{(j)i)sm(l)icos(p2+ 

(2vr)^ Jsn-i Jo ^ (g) 

L(r)M'((/)i)r^^ cos (pi cos </>2) r'^'^drda{$). 



da 



We can then upper-bound these integrals in terms of Fq (the radial component of /). See 
Appendix B.l for details. The final result is: 



E{b^{I -9F)b\a<r]) 

< ^ n-l 
- Pr[a < r?] (27r)2 ° 



/>oo 

i(n-2) / Fo(r)V"-3dr 

Jl/{\rie) 

^ Jl/{Xve) -^l/CAr^e) ^ 

+ (2n + 9 + ^)eA / Foirfr^-^dr 

Jl/(Xve) 
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In sections 4 and 5, we will explain how this is used. 

3.4 The variance of b parallel to 6: In a similar way, we can upper-bound the variance of b ■ 9. 
See Appendix B.2 for details. 

4 The Ball in 

Let be a ball in M", of radius (3. In this section we will analyze the curvelet transform of the 
function f{x) = [1/ y^vol(i?) if x G i?, otherwise]. This is the wavefunction one gets by quantum- 
sampling over B. 

We assume n > 4, and we use the window functions W and V specified in Section 3. We set the 
parameter A to lie in the range 2tt (3e/n < A < 2 • 27r/?e/n. We show the following: 

Theorem 3 Almost all of the power in f is located at frequencies \k\ > 1/A; J^^^i/x < 
For any rj < the probability of observing a fine-scale element a < rj is lower-bounded by: 

Pr[a < rj\ > ^4(1 — Furthermore, if rj < (l/2e^)(l — then the variance of b, in the 

directions orthogonal / parallel to 0, conditioned on a <rj, is upper-bounded by: 

E{b^{I - ee^)b I a < ??) < r?/32(14300 + 0(i)), (8) 

E{{b-ef\a<r,)<p\2A2 + 0{l)). (9) 

The first claim shows that only an inverse-polynomial fraction of the probability mass lies below 
the low-frequency cutoff; this justifies our use of the curvelet transform and Theorems [1] and [2j The 
second claim shows that, for any sufficiently small constant 77, we observe scale a < rj with constant 
probability. This is due to the fact that / has a lot of power at high frequencies (a "heavy tail"), 
which is caused by the discontinuity of / along the surface of the ball. (For comparison, one would 
not observe this behavior if / were, say, a Gaussian.) The third claim shows that, when a < r], b lies 
within distance 0{^/rj(3) of the line that passes through the center of the ball with direction 9. Also, 
b lies within distance 0(/3) of the center (and we expect, though we do not prove, that this distance 
is also lower bounded by 0,(13)). 

As mentioned previously, it is remarkable that these bounds do not depend on the dimension 
n. It is also interesting that this concentration of probability mass can hold even when the window 
function g only C^-smooth. By contrast, in order for to be asymptotically rapidly decaying, 

must usually be C'^ or C°°-smooth. 

We prove this using our results from Section 3. Note that the curvelet transform behaves in a 
simple way when we translate the function /: if g{x) = f{x — z), then using equation ([T|), we see 
that Tg(a, b,6) = T f{a,b — z, 6). Thus, without loss of generality, we can assume that the ball B is 
centered at the origin. In this case / is a radial function. 

We will be interested in the Fourier transform of /. We write f{x) = /o(|2;|), where foir) = [C 
if r < /?, otherwise], and C = Xj \J yo\{B) = (5o/3"/n)~^/^, where Sq is the surface area of the unit 
sphere in M". Then the Fourier transform of / is given by f{k) = FQ{\k\), where 

P^^P) = y/i^Jn/2{27rpP), (10) 

using the definition from Section 3, and the identity ■^{^'^ Ju{z)) = z^ Jy^i{z) (see [T], eqn. (9.1.30)). 

The behavior of Fq{p) depends on the behavior of the Bessel function Ju{z) when 1/ > and 
2; > 0. Ji,{z) is very small when z <C z^, it undergoes a transition near z ^ v, and it is approximately 
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given by ^^2/{-kz) cos{z— ^utt — ^tt) when z ^ u (see [T]). Thus, our intuition is that Fo{p) when 
p < (n/2)/(27r/3), and Fo{p) « (const)/ p^"-'^'^^/'^ times an oscillating factor when p > (n/2)/(27r/3). 

We now sketch the proof, omitting the details. First we choose A ~ (const) ■ (27r/3) /n, so that very 
little power lies at frequencies below 1/A. Substituting into ([2]), we get that Pr[a < r]] (const) ■ rj. 

Then we bound the variance of h perpendicular to 0, as follows. (A similar argument holds for 
the variance of h ■ 6.) We start with ([7]), and use the fact that r > l/(Xr]e) implies r ■ X-qe > 1: 



E(P'(I 

< 



-ee^)b\ a<r]) 
1 n-1 



Pr[a < rj] (27r)2 



i(n-2)(A7?e)^ 



Fo(r)2r"-idr 



Xrie 



l/{Xve) 
oo 



+ (2n + 9 + ^)eA(Ar/e) / 

Jl/i\ve) 



1/iXrie) 
10 



17(Ar?e)3 

^ -l/CAr^e) 

Fo(r)V"-ifir 



A standard calculation shows that -Fo(r) = —y^^ "'^^ <^(n/2)+i(27r/3r)(27r/?), which behaves simi- 
larly to FQ(r), except that it oscillates differently and is larger by a factor of 2ttP. Thus we can get 
a reasonable estimate by replacing FQ(r) with Fo(r)(2iTP) in the integral: 



E(6^(/ 

< 



-ee^)b\ a<7]) 

1 n-l 

Do • 



Pr[a < 7]] (27r)2 



i(n-2)(A7?e)^ 



X 



+ (2n + 9 + ^)eX(Xrje) 



l/(Ar,e) 



Fo(r)2r"-Mr 



Using the definition of A, 

E(^(I 

< 



-9e^)b\a<7]) 
1 n-l 



(const) 



n. I 



Pr[a < rf\ (27r)2 
+ l(corest) ^ ^' rj + (const) ^ ^' rj 1 + (const) ^ ^' rj 



Fo(r)^r'^~'dr. 

l/{Xrie) 

Now recall the expression for Pr[a < rj] given by ([2]). We expect the two integrals to roughly cancel 
out, so we get 

E(t' (I -6e^)b\ a <rj)<0(/3S). (11) 

This argument can be made rigorous, using known results about Bessel functions Ju(z). However, 
there is a technical obstacle: our theorem concerns the case where z is roughly proportional to v. 
This is still in the transition regime, and the usual asymptotic expansions for Ju(z) do not work 
here (they only work when z > u'^, or when z/u is some fixed ratio). Fortunately, there are useful 
bounds on the quantity Mt^(z) := \J Ju(z)'^ + Y^(z)'^, and representations of Ju(z) and Y,y(z) in terms 
of a modulus and phase, that do work in this regime [H [30]. This leads to a rigorous proof of our 
theorem — see the Appendix for details. 
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5 Spherical Shells 



We now consider the curvelet transform of a function supported on a thin spherical sheh in M" . We 
win show results similar to the previous section, except that they now depend on the thickness of 
the shell. Intuitively, when the shell is very thin, we can measure very fine-scale elements (a small) 
with significant probability, and b is tightly concentrated around the wavefront set. 

Without loss of generality, we can assume the shell is centered at the origin (see Section 4). 
So consider the following function on R", f{x) = [C if P < \x\ < P + 5, otherwise], where C = 
-\- 6)"-Bq — P"-Bq, and Bq is the volume of the unit ball in M". This represents a uniform 
superposition over a spherical shell centered at the origin, with inner radius (3 and thickness 5. We 
call this a spherical shell with "square" cross-section. 

This is the exactly the kind of state that appears in our quantum algorithm. However, it is 
difficult to analyze, as its Fourier transform involves a linear combination of two Bessel functions 
oscillating at different rates. We are interested in the case where 6 <^ p. In Appendix D.l, we give a 
heuristic explanation of why the curvelet transform of this state will be tightly concentrated around 
the wavefront set. (This holds when 5 < P/n.) 

Here, we give a more rigorous argument, for spherical shells that have "Gaussian" cross-sections — 
when (5 <C /3, these functions are similar to the above, but they are analytically tractable. We 
define / = Cfg * q, where: Cj is a normalization factor; g \s a. Gaussian of width 5, that is, 
g{x) = (^""/^ exp(— 7r|xp/(^^); q is the measure supported on the sphere of radius jS around the origin, 
which is obtained by restricting the usual volume measure on M"; and the star denotes convolution. 
Intuitively, q represents a shell with infinitesimal thickness, and / represents a "smoothed" shell 
with thickness 5. 

The Fourier transform of / is given by / = Cfg ■ q, where: g is a Gaussian of width 1/5, g{k) = 
5"/2exp(-7r(52|^|2); and q is given by q{k) = Qo{\k\), where Qq{p) = ^^^^ J(n/2)-i(27r/3p)/3'^/2. 
Intuitively, the Fourier transform of the spherical shell is somewhat like the Fourier transform of the 
ball, except that it decays more slowly (i.e., has more power at high frequencies), for frequencies up 
to roughly 1/(5; the power at frequencies above 1/5 is suppressed by g. 

Note that this is quite similar to equation (jl28p , describing a spherical shell with "square" cross- 
section, when we substitute in the upper and lower bounds on C/ (to be given later in this section). 
This suggests that the spherical shell with "square" cross-section can indeed be approximated by 
one with "Gaussian" cross-section, when b P/n. 

We will prove bounds on the continuous curvelet transform, with the same window functions as 
in Section 4. We use a slightly different scaling parameter A: we set A = 7^^, where /3 is an estimate 
of the true radius of the shell, which satisfies [3 < (3 < S(3, for some S > 1. We assume that the 
dimension n is at least 4, and we assume that the thickness of the shell is small compared to the 
radius: 5 = e(5, where e < min^ , (Note that for these values of n, ^ < and 

^ < ^, so the second and fourth conditions are actually redundant.) Under these assumptions, we 
prove the following: 

Theorem 4 Almost all of the power in f is located at frequencies \k\ > 1/X: J^^^i^x — §• 

Let rjc = {5/ f3){n — 2)/e. The probability of observing a fine-scale element a < r]c is lower-bounded 
by: Pr[a < rjd > 0.045. Furthermore, the variance ofb, in the directions orthogonal / parallel to 9, 
conditioned on a < rjc, is upper-bounded by: 

E{b^{I - ee^)h\ a<rjc)<{n- l)e/32(507 + 0(i)) • S, (12) 
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E{{b-ef\a<il,)<^\22, + 0{^)). (13) 

The proof uses a similar strategy to what we showed in Section 4. The intuition is as follows. First 
write f{k) = -Fo(l^l)- The important difference (compared to Section 4) is that here FQ{r) decays 
more slowly, like l/r*^"^"*^^/^, for r < 1/5. Substituting into ([2]), we see that a constant fraction of 
the probability mass lies at frequencies of order 1 /6. So with constant probability, we can observe 
fine-scale elements a < rjc where rjc is of order 6/\. Note that rjc ^ 5/X ^ {const)[5 / f3){n — 2) ~ 
{const)e{n—2). So, when the shell is very thin, r/c will be very small, and b will be tightly concentrated 
around the wavefront set. 

The rigorous proof is given in Appendix D.2. 

6 A Fast Quantum Curvelet Transform 

6.1 The Discrete Curvelet Transform: First, we describe the discrete curvelet transform, which 
has been studied in the classical setting The discrete curvelet transform takes a function f{x) 
and returns a function T f{a,b,9), where both functions are defined over finite domains. This is 
constructed analogously to the continuous curvelet transform, except that one now uses the discrete 
Fourier transform on (Zjvf)", and a discrete set of scale/direction pairs {a, 9). 

The discrete Fourier transform is defined as follows. We assume that / is defined on a domain 
Z consisting of a discrete grid in a finite region of M". For example, let Z = (aZ)" n [— -L,-L)"', 
the intersection of a tightly-spaced square lattice and a large cube. Also let Z = (^Z)" n 
[~2ct' ^ct)"' '^^^ discrete Fourier transform maps / to a function / defined on Z, as follows: 
/(^H'(^)"/'E.-ez/(^)e-'"''-^ and /(f) = (^)"/2 ^.^^^ /(^)e2-*-. 

One can argue that this is approximates the continuous Fourier transform in the following sense. 
Let fcont be a function on M", and let fcont be its continuous Fourier transform. Suppose that fcont 
is supported inside the cube [— L,L)", and fcont has all except an e fraction of its probability mass 
inside the cube [~^) 2^)"' Then there exists a function f^is on Z, with discrete Fourier transform 
fdis, such that fdis ~ cr'^^'^fcontlz and fdis = (^)"'^^/contlz' errors whose total probability 

mass is roughly e. See Appendix E.l for more details. 

Now the discrete curvelet transform is given by 

Fjia, b, 6) := (^)"/2 /(^);^^_^^(fc)e2-^^-^. (14) 
feez 

The "location" variables x and b take values in Z, and the "scale" and "direction" variables (o, 6) 

take values in some discrete set G, which we will describe below. The window functions y ^ are 

defined over Z, and are constructed so that ~ ^ ^- This ensures that the curvelet 

transform can be realized as a unitary operation on the space spanned by the states \k, a, 6). 

Recall from Section 2 that each window function q is supported on a "sector" ^ that has 
angular width -y/a, inner radius 1/a and outer radius 2/a. To satisfy the above condition on Xa 

want to choose a discrete set of values (a, 6*), that corresponds to a discrete collection of sectors 
that forms a "tiling" of the frequency domain. Intuitively, this is done by setting a equal to powers 
of 2, and sampling from a mesh with angular spacing ^/a on the sphere Sn-i- Then the sectors 
fit together nicely, as in Figure [TJ (This picture is a slight oversimplification; actually, since we 
want the window functions X a to decay smoothly to zero, we should make their supports overlap 
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Figure 1: Construction of a 2-D discrete curvelet transform — "tiling" of the frequency domain into 
sectors plus special sectors Siow and Shigh for very low and very high frequencies. Within the 
annulus, a takes on values 1/2'^, for k = 0,1,2,3,4. Each sector has inner radius 2^ = 1/a, outer 
radius 2^+^ = 2/a, and angular width (7r/2)/2L^/2J (7r/2)^. 

slightly.) We will describe a construction of this kind in the next section; other constructions were 
given in [BJ [20] . 

This discretization affects the values of a and 9, relative to the continuous case. Intuitively, the 
"discrete" a can differ from the "continuous" a by a constant factor, and the "discrete" 9 can differ 
from the "continuous" 9 by an additive error of size ^/a. 

6.2 The Quantum Curvelet Transform: The quantum curvelet transform is the unitary opera- 
tion that maps 

^/(x)|x)|0,0) ^^r^(a,6,0"')|6)|a,e"). (15) 

^ a,b,e 

This can be implemented as follows: first apply the quantum Fourier transform (QFT), then the 
operation X that maps |A;)|0, 0) ^ \k) q-Xq g(^)|fli ^)) aiid then the inverse QFT. 

We want to compute this in time polynomial in n and log(M) (where M = 2L/a is the length of 
the discrete Fourier transform). This is possible for the QFT. But it is not clear how to perform the 
operation X, for a generic choice of the window functions Xao- (Note that we want the functions 

fj to be C^-smooth, so their supports will necessarily overlap; thus the operation X must prepare 
a superposition containing 

2e{n) terms.) 

Nonetheless, we can perform X efficiently in two cases: (1) when the window functions are 
indicator functions supported on disjoint sets, and (2) when the window functions are smooth "bump" 
functions that can be expressed as products of 1-D functions using spherical coordinates. The first 
case has poor analytic properties, but the second case is a reasonable approximation of the curvelets 
used in Sections 3-5. Thus we get an efficient quantum curvelet transform. See Appendix E.2 for 
details. 
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7 Quantum Algorithms using the Curvelet Transform 



7.1 Single-shot measurement of a quantum-sample state: Consider the following problem. 
Let S be a ball of (unknown) radius f3 centered at some (unknown) point c in R", for n > 4. We 
are given as input: n, the dimension; /?, an estimate of the radius of the ball (we arc promised that 
/?/2 < P < P)', R, an outer bound on the location of the center (we are promised that |c| < R); fi, 
the desired accuracy of our answer; a description of the set of grid points G = (aZ)" n [— L,L)", 
such that L > R + P and a < i^ggob+QT ^^^^ some constant Qi, to be specified later); and a 

quantum state —7= — = '^g^nnB ^^^^ ^ single quantum-sample over the ball B. We are then 

asked to output a point z in R", that lies within distance /i of the center c. 

We propose the following algorithm. Intuitively, this algorithm uses the curvelet transform to 
find a line that passes near the center of the ball, then guesses a random point along this line. 

Algorithm 1: 

Let IV') be the input quantum state. 

Set ?7 = i4300+\Qi/n) ' where Qi is some constant. 

Apply the fast quantum curvelet transform, with A = Smin = 1) Smax = Ig ^ + 3. 

Measure the scale a = 2~*, location h and direction 0. 
U a > rj, then return "no answer." 

Set = \/3a/242 + {Q2/n), where Q2 is some constant. 
Guess some u e [— 1, 1] uniformly at random. 
Return the point b' = h + uQ^(59. 

We are especially interested in instances where the error /Lt is a constant fraction of the radius (3, 
i.e., II = for some fixed v < 1. Wc conjecture that, for any Algoritlim 1 solves these instances 
with probability independent of tlic dimension n. (In other words, tlic success probability lias 

a "heavy tail." ) This is a sharp contrast to what happens in the classical case: if we choose a single 
point uniformly at random from the ball, then the success probability is z^", which is exponentially 
small in n. This is because, in high dimensions, most of the volume of the ball lies near its surface. 
This is bad for classical sampling, but it helps the quantum curvelet transform, which works by 
finding a line normal to the surface of the ball. 

We now show how our results from Section 4 support this conjecture. We prove the following: 

Theorem 5 Consider a "continuous" analogue of Algorithm 1, using the continuous curvelet trans- 
form ouerR". This algorithm succeeds with constant probability > Jl(z/^). 

We will also argue, non-rigorously, that the discrete algorithm will behave like the continuous one, 
provided that the grid G is sufficiently fine. When the grid is chosen properly, the discrete algorithm 
runs in time poly(n, log i?, log /3, log ^). See Appendix F.l for details. 

We remark that it should be possible to achieve a better success probability, Q.iu'^), using the 
quantum curvelet transform. Here, we showed that b was within distance 0(/3) of the center, so in 
the last step of the algorithm, we simply guessed a point along the line, with success probability 
^{v). But in fact, b should lie at distance ~ (3 from the center, so we should be able to guess one of 
the two points 6± (39, with success probability 0,(1). 

We also remark that classical sampling becomes more powerful if one is allowed any constant 
number of samples, instead of just one. By sampling k random points from the ball and taking 
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their average, one can find the center with accuracy :^[)/^/k, with constant probabihty. (However, 
for fixed k, the success probability does not have a "heavy tail," i.e., one cannot expect to get 
better accuracy with significant probability. This is because, in high dimensions, random sampling 
produces k vectors that are nearly orthogonal.) 

7.2 Quantum algorithm for finding the center of a radial function: Let / be a radial 
function on M" (where n > 4), centered at some point c, and taking values in some arbitrary set. 
Suppose that the level sets of / are concentric spherical shells of thickness 5 centered at c, i.e., / is 
constant on each shell, and / takes on distinct values on different shells. (Note, in previous versions 
of this paper, we made an additional assumption, that one can efficiently compute the radius of a 
shell, given the value of / on that shell. This assumption is no longer needed.) 

Consider the following problem. We are given as input: n, the dimension; i?, an outer bound 
on the location of the center (we are promised that |c| < i?); 6, the thickness of the spherical shells; 
/i, the desired accuracy of our answer; and an oracle that computes the radial function /. We are 
asked to output a point z in M", that lies within distance fi of the center c. 

We propose the following algorithm. The basic idea is to prepare a quantum superposition over a 
large ball around the origin, then measure the value of / to get a superposition over a spherical shell 
centered at c, then apply the curvelet transform, and find a line that passes near c. The algorithm 
does this twice, then returns the point on the first line that lies closest to the second line (note that, 
with high probability, the two lines are nearly orthogonal). 

Algorithm 2: 

Set R' = nR. Let B he a ball of radius R' around 0. 

Choose the grid G = (uTf n L)", where L = 2R' and a = 5/400. 

For i € {1,2}, do the following: 

Prepare the state ^^^^^^ Ylx&GnB 1^)' using the methods of [2] or [T7] . 

Compute the value of / in an auxiliary register, and measure it; call this y^^\ 

Set =R' -R. 

Set r/W = ((5//3W)(n- 2)/e. 

Apply the fast quantum curvelet transform. 

With - „_2 , - •Smax - Ig ^ + ^. 

Measure the scale a^*^ = 2"'^'*'', location 6^*^ and direction 0^^\ 
If a^*) > ry^*-*, then return "no answer." 
End for. 

If |^(^) • ^(^)| > 3/4, then return "no answer." 

Set r = . e(^), s = . (gi) -b(^)),t = e(^^ • {bW _ 6(2)). 

Return the point + b^^l 

We conjecture that Algorithm 2 finds the center with arbitrary precision fi, provided that 5 is 
sufficiently small, i.e., the radial function / computed by the oracle is sufficiently "precise." Let us 
assume that 6 < ■ f^^^^i^ij^ ■ 76i+(Qi/n) ' some constant Qi to be defined later. We conjecture 
that Algorithm 2 then finds a solution with constant probability, independent of the dimension n. 
Thus, only 0(1) oracle queries are needed. This is an improvement over the classical case, where 
i7(nlog^) queries are required. (The indicates that we are omitting some log factors.) 

We now show how our results from Section 5 support this conjecture. We prove the following: 

Theorem 6 Consider a "continuous" analogue of Algorithm 2, using the continuous curvelet trans- 
form overW^. This algorithm succeeds with constant probability. 



14 



We will also argue, non-rigorously, that the discrete algorithm will behave like the continuous one, 
provided that the grid G is sufficiently fine. When the grid is chosen properly, the algorithm runs 
in time poly(n, log ii, log ^). See Appendix F.2 for details. 

7.3 Classical lower bound: We claim that any classical algorithm for finding the center of a 
radial function must use at least r2(nlog ■^) queries. (The indicates that we are omitting some log 
factors.) 

Our intuition is as follows. Any algorithm can be described as a decision tree, where each node 
represents a query to the oracle /, and the algorithm chooses which branch to follow depending on 
the oracle's answer. However, the values of the function / are meaningless by themselves, so when the 
algorithm receives an answer from the oracle, the algorithm cannot do anything besides comparing 
this answer with the answers returned previously. Thus, after its /c'th query, the algorithm can 
choose one of at most k distinct branches. 

It follows that, if the algorithm makes £ queries, the number of possible outputs (i.e., the number 
of leaves in the tree) is at most £!. In order to solve this problem, however, the algorithm must be able 
to output at least (R/fj,)"" different points. So we have {R/ fi)"' < £\, which implies £ > 0,{nlog 

A formal statement and proof of this result is given in Appendix F.3. 

7.4 Finding the center through multiple iterations: We now describe a variant of Algorithm 
2, for finding the center of a radial function. This new algorithm will use multiple iterations, and a 
larger number of queries, but it has a less demanding requirement on the thickness of the shells that 
form the level sets of the radial function /. 

First, we describe a single iteration of the new algorithm. We call this procedure OneRound(). 
This is similar to Algorithm 2, but it starts out with a promise that the center lies within distance R 
of some point p, and it returns a point q that lies within distance ^/R-\/ [ijl of the center. OneRound() 
also takes a parameter 5 > 1 that controls the accuracy and success probability: OneRound() returns 
a point q (instead of "no answer") with constant probability, and when this happens, the point g is 
accurate with probability > 1 — O(^). 

Procedure OneRound(i2, p, S\. 

Set R! = nSR. Let i? be a ball of radius R' around 0. 

Choose the grid G = (ctZ)" n [-L, L)", where L = 2R' and a = (5/400. 

For i € {1,2}, do the following: 

Prepare the state ^^^^^^ X^^eGnB 1^)' ^sing the methods of [2] or [T7]. 

Define the function g{x) = f{x + p). 

Compute the value of g in an auxiliary register, and measure it; call this y^^\ 

Set =R' -R. 

Set r/» = {6/i5^^){n-2)/e. 

Apply the fast quantum curvelet transform. 

With - „_2 , Srrtm " ^^nax - Ig ^ + 3. 

Measure the scale a^*) = 2~**''\ location 6^*) and direction 0^'^\ 
If a^*-* > T]^*-*, then return "no answer." 
End for. 

If 1^(1) • ^(2) I > 3/4^ then return no answer. 

Set r = • e^^), s = • (6(1) - ^W), t = 9^^) ■ (gi) - m). 

Return the point q = p+ =j^9'^^'^ + b^^'^ . 
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Now we describe the full algorithm, with multiple iterations. This algorithm begins with a point 
at distance R from the center, then uses OneRound() to find a point at distance ■\/R^/J^/2 from the 
center, and by repeating the procedure, shrinks the distance to i?^/^(/i/2)^/^, R}/^{ii/2y/^ and so 
on. It may seem surprising that the distance decreases by more than a constant factor during each 
iteration. Intuitively, this is because the spherical shells used by the algorithm arc not exact dilations 
of each other. Recall that the shells have different radii /3, but they all have the same thickness 5. 
The larger the radius /?, the smaller the ratio e = 5/13] so a larger shell allows a significantly more 
precise determination of its center. In a sense, the algorithm makes more progress during the early 
iterations, when the spherical shells are larger. 

Algorithm 3: 

Set Rcur = R and Pcur = 0. 

Set Hiter = [Iglg^l, S = {9A)niter and ntries = 910 log S. 
While Rcur > 1^ do: 

Try running OneRound{Rciir, Pcur, S) up to ntries times. 

If OneRound() returns "no answer" on every attempt, 

then return "no answer." 

Let q be the point returned by OneRound() on one of the successful attempts. 

Set Rcur = VRcury/JjJ^ and Pcur = 
End while. 
Return Pcur- 

We conjecture that Algorithm 3 will succeed when 



128(10 lglg^)3n2(507 + ^)' 

which is a weaker requirement than that of Algorithm 2, where 5 had to scale like 1/R. We conjecture 
that this algorithm then finds a solution with constant probability. Note that this algorithm uses 
0(lglg ^ Iglglg ^) queries, which still beats the classical lower bound of r2(n log ^) queries. 
We now show how our results from Section 5 support this conjecture. We prove the following: 

Theorem 7 Consider a "continuous" analogue of Algorithm 3, using the continuous curvelet trans- 
form over R". This algorithm succeeds with constant probability. 

We will also argue, non-rigorously, that the discrete algorithm will behave like the continuous one, 
provided that the grid G is sufficiently fine. When the grid is chosen properly, the algorithm runs 
in time poly(n, log R, log ^). See Appendix F.4 for details. 

8 Conclusions 

We introduced the curvelet transform as a tool for quantum algorithms, and demonstrated how it 
can solve problems involving geometric objects in M". We showed that: (1) for functions with radial 
symmetry, the contim.ioT.is curvelet transform concentrates probability mass near the wavefront set; 
(2) a quantum curvelet transform (which is a discrete approximation of the continuous curvelet 
transform) can be implemented efficiently; (3) this leads to quantum algorithms for approximately 
finding the center of a ball in M", given a single quantum-sample state, and for exactly finding the 
center of a radial function in R", using 0(1) oracle queries. 
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There are several ways in which these results might be extended. Perhaps one can adapt these 
quantum algorithms to solve more general problems, like finding the center of an ellipsoid. Perhaps 
the quantum speed-up can be amplified using a recursive construction, as in [Hill]. 

A general open problem is to understand the behavior of the curvelet transform on more com- 
plicated shapes. Can one prove that the probability mass of the curvelet transform is concentrated 
near the wavefront set, for arbitrary functions on M"? That would generalize the results of this 
paper, [lOj and [7]. Also, can one rigorously bound the approximation of the continuous curvelet 
transform by a discrete one? 

Another problem is to find new quantum algorithms based on the curvelet transform. For ex- 
ample, can one construct a curvelet transform over F^, that could help to solve hidden polynomial 
problems [12j? Are there quantum states with "wavefront" features, from which the curvelet trans- 
form could extract useful information? Some candidates are quantum-sample states over convex 
polytopes [2, 17J, and states produced by the evolution of a quantum walk. 

One might also try to use the output of the curvelet transform in a more sophisticated way. In 
this paper, we simply measured the output state, and we made very little use of the scale variable 
a, which measures the "sharpness" of the wavefront discontinuity. 
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A The Continuous Curvelet Transform 

We now sketch the proofs of Theorems [T] and [21 
First, for any a and 9, let us define 

9a,e<^)'-= (X,br fKAe<^)db. (17) 

We claim that 

9a/k) = \x,/k)\'nk). (18) 

To see this, write 

(w/) = X,j:o/(2;-6)/(yidy^= i%,o,e*fm, 

where we define 

Thus we can write g{x) as 

Taking the Fourier transform, we get 

^a,.^^) = 7,,o>-(^)(7:,o>-)"(^)/(^) = KoM'hh 

using the fact that 

This proves equation ([18]). 

Next, we claim that, for all |A;| > 1/A, 

r/ = 1- (19) 

To see this, proceed as follows. Fix A;, and write 6 in spherical coordinates centered around k, 
such that 6i is the angle between 6 and k. Then we have 
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Then substitute into the integral: 

da 



^5"-i 

da 



JS"-i ^ sm(t'i) / ( 



w{\a\k\fv{ei/./^fs\n^-'^{ei/^/^)deid(T{e2,...,en-{ 



3/2 ' 



« 

Note that V{9l/^/a) is nonzero only when 9i £ [0, vr-^/o], so we can restrict the integral to this range. 
Then change variables, 9[ = 9\l\fa^ to get: 

J5"-i 



C [ W{Xa\k\fV{9[f sm''-\9[)d9[V^ da{92, . . . ,9n-i] 
Jo Js"-^ 



da 

3/2 



a 



W{Xa\k\fV{9ifda{9) — 


0^ Js"-^ 



1, 



using the admissibility conditions. This proves equation p9 
We now prove Theorem [TJ We write 



Tfia, b, 9)-f^^ff{x)dfi{a, b, 9) 

a<l ' ' 

, / i^a,br fhaAe<^)dbda{9)—'' 

1 Imn a 



JS 

/ 9a/^)dai9) 

JS"-! ^ 



+1 

da 



and we claim that this equals /(x). Taking the Fourier transform, we get 

and we claim that this equals f{k). Using equations (llSp and ()19p . we rewrite this integral as: 

i2 S/^N . 



is"-i « 

(The last equality holds because of (fT9]) when |A;| > 1/A, and because f{k) = when < 1/A. 
This proves Theorem [TJ 

Finally, we prove Theorem [2j We write 

Tf{a,b,9)\^dfi{a,b,9)= [' [ [ ^^J)\^dbda{9)^. 
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We rewrite the innermost integral, applying some of the identities used to prove (jlSp : 

Substituting in, and using (fT9|) . we get: 

|r/(a,6, e)\'^dfi{a,b,e) 



a<l 

1 



|/(A:)|^dfe= / \fix)\^dx. 



J5 
2 



This proves Theorem [2l 

B The Curvelet Transform of a Radial Function 

B.l The variance of b perpendicular to 9 

This is a continuation of Section 3.3. Recall that we have L{r) = Fo{r)W{Xar), and M{(j)i) 
^(<Ai/\/fl)^a(0i)- Then we write: 

\^^^^^^{k)\^dk 

L'{r)M{(f>i) sm(f>i cos (j)2+ 



Jo 

N 2 

L{r)M' {cf)i)r^ cos 01 cos 02 ) drda{(j)) 



(20) 



= lArlAlh + '^ISrlBlh + Icrlcih, 

where we define 

lAr= L'irfr'^-^dr, Iai= M(0i)2 sin'^ 0i#i, (21) 



JQ 



lBr= L'{r)L{ry-^dr, Ibi = M'(0i)M(0i) cos 0i sin"" ^id^i, (22) 
Jo Jo 



CO 

Icr= I Hrfr^'-^dr, Ici = I M'(0i)2 cos^ 01 sin"-^ (/)i(i0i, (23) 
'0 Jo 



'2 



COS^02dcr(02,---,0n-l)- (24) 



Sn-2 
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This shows that the variance of b perpendicular to 9 is: 



So [ (lArlAlh + '^ISrlBlh + Icrlcih^ 



(2vr)- JO ^ 

A similar formula gives the variance conditioned on observing a <r]: 

1 



a 



n+l ' 



E{iF{i -ee' )b\a < 7]) 



Fr[a < r]] (27r)2 



fi — I rv X \ da 

So / ( lAr^Alh + '^iBrlBlh + Icrlcih 



n+l ' 



(25) 



(26) 



We will now upper-bound the various integrals appearing on the right hand side of equation 
B.1.1 

We begin with the integral lBrlBil2dO'/ci'^^^ ■ A straightforward calculation shows that 

1 



n — 1 



da {4)2 



(27) 



We can simplify Ibi, by integrating by parts: Ibi = — Jq ^M{(j)i)'^{—1 -|-ncos^ 4>i) sin"~^ ^ic^^i- 
Substituting in the definition of M{(/)i), changing variables, and using the fact that V is supported 
on the interval [0,7r/2], we get: 



Ibi 



V smcpi J 

iy((/>i/y^)2a(2"-i)/2 sin(</.i/V^)"-2(-l + ncos^ 0i)#i 



/ iy(LJi)^a"'sin""^cJi(-l + n cos^(wi V«))c?^i 
Jo 



tt/2 



iy(u;i)^a" sin" ^ciJi(l — n cos^(LL'iV^))dti;i, 



hence 



I-k/2 

\Ibi\ < ia"(n- 1) / y(wi)2sin"-2^ida;i 
Jo 

Combining this with I2, we get: 



\lBih\ < ^a" / V{L0i)^da{uJucP2, • • • , K-i) = ^a"- 



(28) 



Next, we can simplify iBr, by integrating by parts: iBr = — Jq" ^L(r)^(n — 2)r" ^dr. Combining 
this with Ibi and ^2, substituting in the definition of L(r), and exchanging the integrals, we get: 



^BrlBlh-^ 



(r)2(n-2)r"-3dri 



JO « JO Jo 

n f'°° dn 
= / / iFo(r)2VF(Aar)2(n - 2)r"-3dri — 
Jo Jo « 

r°° n dn 

< / ^Fo{rf / Ty(Aar)21 — (n-2)r"-3rfr. 
Jo Jo 0, 



da 
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By the definition of W, 



and vanishes when r < l/{\rie). Thus we have: 



rv 
Jo 



da 



,n+l 



< 



lin-2) [ 
Ji 



Fo(r)V"-3dr. 



l/(A»?e) 



(29) 



(30) 



B.1.2 



Next, consider the integral Jq lArlAihda/a^^^ . We already have a bound for I2. For 7^1, we write: 

lA. = /V(<^,/v^)^a(--)/^(^^^^i^)"-%in"<^,#, 
Jo ^ sin(/)i / 

= / F(u;i)^a"'sin"~^wisin^(u;ii/a)dwi. 
Jo 



Using the fact that V is supported on [0, 7r/2], and the simple bound sin^(a;iv^) < Lofa < \a, we 
get: 



p-n/2 

Jo 



Combining with I2, we get: 



< lAlh < / V{uJifda{uu(t)2. <t)n-l] 

4 n - 1 Jsn--L 

We now turn to Iat- First, combining with Iai and 12, we have 



7r2 

4 n- 1' 



< / lArlAlh-^^ < I lArda 



rv 
Jo 



rv 
I 



4(n - 1) 



(31) 



(32) 



We can upper-bound lAr as follows. Note that, for any two functions, the Cauchy-Schwarz 
inequality implies that 



+ 2{f,g) + \\gf < Wff + 2||/||||5|| + H^H' < 2\\ff + 2\\gf; 



(33) 



in the last step we used the arithmetic-geometric mean inequality, y/ab < for a,b > 0, with 



and b = \\g\\ ■ Thus we can write 



where we define 



0<lAr = J (^F^{r)W{Xar) + Fo(r)VF'(Aar)Aa)V"-^dr 
= IIG1 + G2IP < 2||Gi||2 + 2||G2||2, 



Gi(r) = F^(r)VF(Aar)r("-i)/2 
G2{r) = Fo(r)VF'(Aar)Aar("-i)/2_ 



(34) 



(35) 
(36) 



23 



Thus we have 



rv 



0< / lArUlh 



da 



< 



TT 



2(n - 1) ■ 



We then want to upper-bound the integrals /(J'llGi pda and /(^||G2|pda. 
For the first one, we have: 



prj pri poo 

/ \\Gifda= / / F^(r)2t^(Aar)V-Wa 
Jo Jo Jo 

poo pTj 

= / ^o(^)^ / W{\arfdar'^-'^dr. 
Jo Jo 



Using the fact that W is supported on [1/e, 1], we can write 

/ W(Xarfda= W(a?da— < W(af — —<—, 
Jo Jo Jo « 

and vanishes when r < l/{\rie). Hence, 



rv 1 f oo 

/ \\Gifda<- / F^(r)V-2dr. 



For the second integral, we have: 

pr) pT] poo 

/ \\G2fda = / / Fo{rfW'{XarfX^a\''-^drda 
Jo Jo Jo 

POO rf] 

= / Fo^rf / W'i\arfX^a^dar''-^dr. 
Jo Jo 

Note that the derivative of W is given by 

^/(^) = 1^"" sin(27r log r)7r/r, 1/e < r < 1, 
lO, otherwise, 

where = yW^- ^o we can write 

rv r^vr 1 1 

/ W'iXarfX^a^da= / l^'(a)2a2da-^ < / CUT^/afa'^da—. = lTr^{l 
Jo Jo X^f J 1/e Xr 

and vanishes when r < l/{Xr]e). Hence, 



rv 1 

/ ||G2f da < Itt^I - D- / Fo(r)V-^dr. 

Jo ^ •^l/(A??e) 
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B.1.3 

Finally, we consider the integral Icrlcihda/a'^'^^ . We already have a bound for I2. For Ici we 
can write: 



0</ci= /"(F'(0i/V^)(l/V^)Aa((Ai) + F(0i/V^)A;(0i))%os2</>isin"-2,/,i#i 

-/n V / (43) 



|C/i + C/2|r <2||?7ir + 2||C/2r, 



where we define 



UiicPi) = y'((/>i/^/^)(l/^)A„(</.i)cos(/.isin("-2)/2</,i (44) 

^^2(01) = y(0i/y^)A',(</'i)cos(/>isin("-2)/2<^i. (45) 

Then 

< Icih < 2\\Uifl2 + 2\\U2fl2- (46) 

We now evaluate ||C/i|p and ||J72|p. 
For ||C/i|p, we can write 

\\U,f = r V'ic^./V^fa-'a^'^-'y' r'^^'^'//^^ ) cos^ </>i sin"'^ <p,dcp, 
Jo ^ sm(</>i) / 

= y'(0i/V^)2a(2"-3)/2 sin"-2((/.i/V^) cos2 
Jo 

= r V'(wi)2a"-isin"-2wicos2(a;i^)dc<Ji 
JO 

fn/2 

< a"-i / y'(u;i)2sin"-2^ida;i. 
Jo 

The derivative of is given by 

y>^^-^^i-2C,cos{t)sm{t), 0<t<7r/2, ^^^^ 
1 0, otherwise, 



where Cy = y Using these formulas and ([T], eqn. 4.3.127), a straightforward calculation 

shows that 

/ y'(a;i)2sin"-2^idc^i = — 1(1--) / sin^-^widwi. (48) 
Jo 36o V n/ Jo 

So we have 

\\U,f<a--'-^(l--) r sin-' u,du;,. 
35o V n/ Jo 

Combining with I2, we get 

< \\Ui fh < a"-i . (1 - i) . [ da{u;uct>2, • • • , <l^n-i) = |a"-^ (49) 

Next W6 GvalucitG ||^2|P* The derivative of A^, is given by: 

A' = a(2«-l)/4/n-2^ pin(0l/^) \ (n-4)/2 / cos((?!>i/ ^g) _ sin((/)i/^) cos((/>i) \ 
° 2 V sin((/)i) / V Y/asin(0i) s'm''{(j)i) / 
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hence 

' cos{4>i / y/a) s\\i{(j)i / ^/a) cos{<j)i) 



and 



sin(_(/)ij 



cos (pi, 



2 



2ir=/ F('^l/V^)V2""l)/2(«^)2sin"-4(<^l/^), ■ .■, ■ I COS-,pia</)i 

V va sm((pi) / 

y cJi V V sm" ^ ^-^^ cos^ wi^/^ dwi. 

Recall that V is supported on [0,7r/2]. For wi in this range, we have the following crude bound: 
(using [1], eqn. 4.3.81) 

cos(wi) sin(cji) cos(6<Jiy^) ^ 1 ^ sin(cji) sin(ci;iy^) ^1^1 



cos(0i/ya) sin(0i/ya) cos((?;<i)\ 2 ^ 



/a sin(i^i-ya) ^/a sin(u;iya) 00l^/a yfa ^/a 

Also, we have cos^(u;iy^) < 1. Hence, 

f-7r/2 



2f< f V{uifa^{^fsin^-\ui){^fdui 

fV2 

= Aa^-\-^f / y(u;i)2sin"-^a;idu;i. 
Jo 



Note that 

r.7r/2 c- /.7r/2 



r V{ujif sin"-4 L^idwi = f 1 + ) /"^ sin"-2 ^^^^^ 

Jo V n-S/Jo 



using the definition of V , and integration by parts. 
So we have: 

r"7r/2 



2r^ 



<4a"-i(n_2)2(i + _5_) r y(cui)2sin"-2u;idL^i. 

JO 

Combining with I2, we get: 



n da ^ \ n 

0< /c./ci/2^<(| + 2(n-2)(l + ^)) Ic. 

We can bound the integral on the right hand side as follows. 

ricr^= r r Fo{rfW{\arfr-^dr^ 
Jo « Jo Jo 



Foirf I W{Xarf^r'^-^dr. 



(51) 



(52) 



< \\U2ri2 < 4a"-^(^)^(l + ^)^ / Viu;irdaiu;u<P2, <Pn-i) 

<a"-i(n-2)(l + ^). 
So, by substituting into ([36|) . we have 

< Icih < 2 • la^-' + 2 • a^-\n - 2)(1 + ^). (53) 
Finally, we turn to Icr- Combining it with Iqi and I2, we have 

da 
« 
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Using the fact that W is supported on [1/e, 1] 

,da f^^i' 



/ w{\arf% = / W{af^\r < / W{af—e\r < eAr, 
Jo Jo " Jo a 

and vanishes when r < l/^X-qe). Hence 

rv An r°° 
/ Icr^ < eX / Fo(r)V"-2dr. 

Jo a JinXne) 



(54) 



dn / \ f°° 

< / Icrlcih^ < ( 1 + 2(n - 2)(1 + ^))eA / Fo(r)V"-2dr. 

JO a ^ \ / Jl/(\r,e) 



(55) 



B.2 The variance of b parallel to 6 

Finally, we seek to bound the variance of 6, in the direction parallel to 6. The analysis is similar to 
the previous case (i.e., the variance of h orthogonal to 9). 
The variance of b parallel to 9 is: 



E((6-^)^ 



JR" 



(h.9f\Tf{a,bMdb^da{9). 



(56) 



We can simplify this by taking advantage of rotational symmetry. Fix a vector u = (1,0,..., 0), and 
for each 9, let i? be a rotation that maps 9 to u. Then 



{R{b)-uf\Tf{a,R{b),u)\^db^da{9). 



(57) 



Then change variables b ^ R ^{b). The integrand is now independent of 9, so we can do the 
integral. We get: 



-> da 



¥.{{b-9f) = So / / bl\Tf{aAu)Ydb- 
Jo Jm" " 

We now introduce some new notation, 

'^>^^4):=Tf{a,b,9), 



.+1 ■ 



(58) 



(59) 



to emphasize that we view this as a function of b. By equation ([T]) , the Fourier transform of ^ is 
given by 

^a,e('^)=f(^)Xa/k). (60) 

And we have, by Plancherel's theorem: 



E{{b-9f) = So r t ^^2,^ ,rM2X da 
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1 d 



2-111 dki 



(61) 



,71+1 ' 



Now, using spherical coordinates A; = (r, 0i, . . . , 0„_i), we write ^a,u{k) as a product of a radial 
part and an angular part: 

^aA^)=L{r)M{(t>i), (62) 
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where 

Then we have 
where 



L{t) = F„{r)W(Xar), M{,fi) = V{^/^AJ^i). 
A4„.(S) = i.„M(*0|^ + LMM'(^,)§|. 



dr 



dcpi sin 



Now we can expand out the following integral: (note that ^a,u{k) is real) 

/ -irKaik) ^dk = [ r (L'{r)M{<p,) cos </>! 
Jmri oKi Js"-'^ Jo ^ 

-L(r)M'(0i)r-^sin(/)i) r'^'^drdai^) 



KArKAlK2 - 2KBrKBlK2 + KcrKciK2, 



where we define 



Kav 

Kai 



)i )" cos" ^1 sin'" " (pid4>i 



poo 

Jo 

^0 

K2= da{(j)2, . . .,4>n-l) 

JS'^-^ 

POO 

/ L\r)L{r)r''-^dr 
Jo 

f'Tl 

I M'((^i)M((/)i)cos(^isin"-Vi#i 
Jo 

roo 

/ Lirfr^'-^dr 
Jo 

M'((^i)2sin"(/)id(^i. 



KBr 

Kbi 

Kcr 

Kci 

Thus we can write the variance of h parallel to as: 

. So 



E((5 



KArKAlK2 - 2KBrKBlK2 + KcrKciK2 



da 

-. n+l ■ 



' (27r)2 7o 

A similar formula gives the variance conditioned on observing a < rj: 
Ei{b ■ ef \ a < ri) 

= , r^\o r(KArKAlK2 - 2KBrKBlK2 + KcrKciK2) 

Pr[a < r)] (27r)2 Jq \ J 
We would then like to bound the various integrals appearing on the right hand side. 



28 



B.2.1 



We begin with the integral JJ' KBrKBiK2da/a"-^^ . 

Note that K2 = (n — 1)I2, while Kbi = Ibi and Kbt = Ibv- Using the argument from the 
previous section, we have: 



da 

KBrKBlK2—^ 



<i(n-l)(n-2) / Fo(r)V 



(76) 



B.2.2 



Next, consider the integral K ArK AiK2da / aP''^^ . We already have a bound for K2. For i^Ai, we 
write: 

Kai = r Vi(t>i/V^fAaih? cos2 01 sin"-2 ^^^^^ 
^0 



= / y(a;i)^a"sin"'~^a;iCOS^(a;i-v/a)(ia;i. 



Using the fact that V is supported on [0, 7r/2], and the simple bound cos^(a;iv^) < 1, we get: 

1-^/2 

0<Kai< Viuifa^'sm'^-^LJiduoi. 
Jo 

Combining with K2, we get: 

< KA1K2 <a'' I V{iOifda{iOi, 02, ... , 0n-i) = a"". 

We now turn to Kat- First, combining with Kai and K2, we have 

rv 
Jo 

Note that i^T^r = ^Ar, so we can upper-bound Kat as in the previous section: 

0<KAr-<2||Gif + 2||G2f, 

where we define 

Gi(r) = F^(r)VF(Aar)r("-i)/2 
G2{r) = Fo(r)VF'(Aar)Aar("-^)/2_ 

Thus we have 



KArKAlK24^T< r KAr-. 



0< / KArKAlK2—^<2 

« ^^0 



^Jo « Jo 



(77) 

(78) 

(79) 

(80) 
(81) 

(82) 



We then want to upper-bound the integrals /J'||Gi Ipc/a/a and /J'||G2|p(ia/a. 
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For the first one, we have: 



da 
a 



n°° rin 
F^{rfW{Xarfr^-^dr— 

I F^(r)2 / WiXarf—r^'-^dr. 
Jo « 



Using equation (f29]) . we get: 



Jo C' Jl/{\-qe) 



{rfr'^-^dr. 



For the second integral, we have: 



CL Jo Jo 



^^^^"^ar?\^a\^-^dr— 



Fo{rf / W'{XarYX^adar^-^dr. 
Jo 



Recall that the derivative of W is given by equation ()40p . So we can write 



W'{XarfX^ada = / W'iafadaA^ < I Cl{^/afada— 



ff) i-Xrjr 

W'{Xar)^X'^ada = 
10 Jo 

and vanishes when r < \/{Xr]e). Hence, 



1 

l/e 



Jo 



B.2.3 



(83) 



(84) 



(85) 



Finally, we consider the integral KcrKciK2da/a'^^^ . We already have a bound for K2. For Kci 
we can write: 



0<Kci= r(y'(0i/V^)(l/^)Aa((/'i) + F((/<i/^/^)A;((/>i))%in"(/.ici,/)i 



(86) 



where we define 



\\Ui + U2f<2\\Uif + 2\\U2f, 

Ui{(^i) = i^'((/>i/V^)(l/V^)A„((/>i)sin'^/2,/)i 
U2{^i) = y((/'i/y^)A;(</.i)sin"/2 (t>i. 



Then 



0<KciK2<2\\U^fK2 + 2\\U2fK2. 
We now evaluate ||c7i|p and ||C/2p. 



(87) 



(89) 
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For ||C/i|| , we can write 



V sm((pij / 
" y (</.i/^/^)2a(2"-3)/2 sin"-2(^^/^) g-j^2 ^^^^^ 

Jo 

<^a'^/ y'(a;i)2sin"-2^i(ia;i. 
4 Jo 

(In the last step we used the bound sm'^{LOl^/a) < ujfa < ^a.) 
Then, using equation (08]), we have 

Combining with K2, we get 

< \\U1fK2 <^a^-^(l--) - [ daiuuH, <Pn-i) = ^a"(n - 1). (90) 
Next we evaluate ||C/2|P- The derivative of is given by equation ([501) . hence 

mi) = F(0i/^/^)a(^-^)/^(^)sin("-^)/^(^,/V^)( ^°^('^^/^) - -^i<^^IVa)--<^^) ^ ^^^^^ 

\ ya sin((pij / 

and 



211^=/ V{c^r/V-afa^'-'y\^f sin"-^(<Ai/^/^) f ^"^^'^^/^^ - ^^"^'^^/v^) ' sin^ 

V -y/a sin((pi) / 

""^^ ^rf \2 nin-2\2 ■ n-4/ n/'COs(wi) sin(cJi) COs(a;i^) \ 2 .3/ ^n, 

y(cJi)V(V) sm" ^(lji) ^ ^—^ sin^(wiV^)c?u;i. 

V va sin(ti;ivaj / 

Recall that V is supported on [0,7r/2]. For loi in this range, we have the following crude bound: 
(using [1], eqn. 4.3.81) 

cos(wi) sin(LL'i) cos(a;i v^) ^ 1 ^ sin(LJi) sin(wiy^) ^1^1 



^/a sin(ix'i\/a) y/a sin(ci;i-y/a) uJi^/a ^fa ^fa 

Also, we have sin^ (wi -y/a) < uj\a < ^a. Hence, 

\\U2f < r^V(c.i)2a"(^)^in"-^(u;i)(-|)2^adu;i 



TT a [—2- 



f7r/2 



^)2 r y(u;i)2sin'^-^cui(ia;i. 
^0 

By equation (I5ip . we have: 

|2 < ^2^n(n_2)2(-L ^ T^' ^(^^)2 sin"-2^^d^^. 

JO 
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Combining with K2, we get: 



0< ||[/2fi^2 <vrV(^)2(l + ^) / V{voifda{u:i,c^2,---An-i) 







Using equation (j29p . we get 

rlri 

/ Kcr- < / Fo{rfr--'dr. 

Jo 0, Jl/(Xne) 



< / KcrKciK^^ < (2^(„ _ 1) + _ 2)2(1 + / Fo{r)\--'dr. 

'0 0- ^ ^ Jl/{Xrie) 



(91) 



<^a"(n- 2)2(1 + ^). 
So, by substituting into ([89|) . we have 

< KC1K2 < 2 • ^a"(n - 1) + 2 • ^a"(n - 2)^(1 + ^). (92) 
Finally, we turn to Kcr- Combining it with Kci and K2, we have 

< I" KcrKc^K2^ < (¥(n - 1) + ^ (n - 2)^(1 + ^^)) [ Kc/^- 

We can bound the integral on the right hand side as follows. 

1"^ rlri n r°° rlri 

/ Kcr- = / / Fo{rfW{Xarfr^~'dr- 
Jo a Jo Jo « 

Fo{rf / W{\arf—r^-''dr. 



(93) 



C The Ball in W 

We prove Theorem [3l 

C.l The low- frequency components 

First, we claim that almost all the power in / is located at frequencies above some threshold 1/A. 
This justifies our use of the curvelet transform, and theorems [1] and [21 for an appropriate choice of 
the parameter A. 

We start by proving an upper-bound on the integral t^^ Ju{t)'^dt, for > 0. Note that (pQ, 
eqn. (9.1.62)) 

|J.(t)|<^ (i^>-i, t>0). (94) 

Also ([T], eqn. (6.1.38)), 

u] > V2^i^^''+he~'' [u > 0). (95) 



Hence 

I 1 ^ _ 
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and 



/ r^JJtfdt< / ( — ) dt = — ( — ) — 



2v 



47ri/2V2z.; • ^^^^ 



This upper bound is useful when z < 2v je. 

We can now calculate the amount of power contained in the low- frequency components of /: 

J\k\<z Jo -^OP" 



n 



2itI3z 



< n 



t-\j^/,{tf . dt (98) 
1 f2-nl3ez\'n 



TTTT? \ n 



Setting z = n/(27r/3e), we get 



\f{k)\'^dk<—. (99) 

|<n/{27r/3e) 



Recall that we set the parameter A so that A > 2tx fie jn. So the region < 1/A} contains at most 
a l/(7rn) fraction of the total power. 

C.2 The decay of J^{x) 

We now prove some technical lemmas on the decay of Jy{x) for x > 2v^ v > 1/2. These follow 
from classical results on Bessel functions [H [30] , though some care is required near the transition 
region at x ~ z^. In particular, the usual asymptotic expansions for J,y{x) only work when x > u'^, 
or when x = av for some fixed constant a. For our purposes, we use an asymptotic expansion of 
Ju{^)'^ + ^uix)'^, that behaves well when x > v. 

C.2.1 

We start by quoting the following result from (f30], p. 447). Define 



M^(x) = Vj^(x)2 + y4x)2. (100) 

Then for all x > > 1/2, 



irx/x^ — u 



This immediately implies an upper bound on Jy{x)'^ , for all x > 2v, v > 1/2: 

M^f < M,{xf < J ^ — • (102) 

TtVx^ - 1/2 TTX ^yS 

C.2.2 

We next prove a lower bound on | J;y(x)|, for x within certain intervals. Note that Ju{x) is large at a 
zero of Yu{x). We will show that (1) the zeroes of Yy{x) are not too far apart, and (2) Ju{x) is large 
in a neighborhood around each zero of Y^{x). 
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To see this, note that Juix) and Yy{x) can be written in terms of a modulus and phase, 

J^{x) = M„{x) cos e^{x), (103) 
Y^{x) = M^{x)sme^{x), (104) 

where M^{x) is as defined above, and 6u{x) satisfies the equation 

(see [Tj, eqn. 9.2.21, and [3D], p. 514). This imphes lower and upper bounds on 6l{x), for all x > 2z/: 

^ < el{x) < 1. (106) 

First, we claim that for any t > 2u, the interval [t,t + ^] contains a zero of Y^{x). To see this, 
write the following, for any 5 >0: 

rx+5 



Ouix + 6)= e^{x) + f 9Uy)dy > e^{x) + ^5. 

J X 



So 9iy{t + -^) > 9,y{t) + -IT. So Oiy{x) must equal an integer multiple of vr for some a; G [t,t + and 

Yi^{x) must vanish at that point. This proves our first claim. 

Second, let be a zero of Y^{x), satisfying (p > 2v. We claim that, for any 5 E [— 7r/2, 7r/2], 

\Jy{(t) + 5)\>My{^ + 5)cos5. (107) 

To see this, write 

I j^(0 + 5)1 = M^{<t) + 5)\ cos e^{(t) + 5)\ = M^{^ + ,5)|cos \e^{ct> + 5)- e^{<p)\\. 

(The last step follows because Oy{(f)) is an integer multiple of vr.) Then note that 

\e,{<p + 5) - = / 9'Mdy < 



Jcl> 

Hence we have 

|cos \Oy{(l) + 5) — ^1/(0)1 1 > |cos \5\ I = cos 5. 
This proves our second claim. 

C.2.3 

Finally, we prove the following lower bound on a sum of squares of Bessel functions: 

Let vi, . . . S [1/2, i^max]- Let t > 2fniax- Then there exists some t' G [t — ^, f + 

27r I 7r 1 



^ + f ] such that 

m 

7t 



^•^-(0'^^- (108) 



k=l 
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Proof: Essentially, we will show that there must exist a point t' where a constant fraction of the 
functions Jy^{t'Y {k = 1, . . . ,m) are large simultaneously. 

Define the interval / = [t,t + ^]. For each /c = 1, . . . , m, / contains a zero of Y,yi,{x), call it cpk- 
Now define the function 



Note that 



COS^(x - (pk) a 4>k - ^ < X < 4>k + ^, 

otherwise. 



2 

Ju^xf > M^^{xfxk{x) > —Xk{x). 

TTX 



Furthermore, define the function 



m 

..Jx) 

k=l 

and note that 

2 



TTX 

k=l 

Define the interval /' = — f , i + + f ] ; this interval contains the support of all of the functions 
Xk{x) {k = 1, . . . , m). Then write 

/■0fc + (7r/2) ^ 



/ u{x)dx = 2. / Xk{x)dx = m 



So there must exist a point t' G /' such that 



1 vr vr 



u{t') > jjj-^ J^^ u{x)dx > - • m • - = — m, 



and the claim follows. 



C.3 The probability of observing a scale a 

Next we claim that / has a heavy tail. This implies that we will observe fine-scale elements (a small) 
with significant probability. 

Again, we start by proving a lower bound on J^°° J^{t)'^dt, when v is of the form m or m+{l/2) 
(where m is an integer), u > 1, and z >2i'. We will show that: 

j t-^Utfit>(l--)^^-^^. (109) 

First, consider the case of = m. We assume m > 1 and z > 2m. Using (pLj, eqn. 11.3.36), we 
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can write 

roo 



^ m— 1 

- [Mtf + Jm{tf + 2 Mtf' 



m—1 

2m 

k=l 
A:=l 

^ m—1 
> —(-Jmizf + 2Y,Jk{zf 
^ k=l 

Then, using the lemma from the previous section, we get the foUowing, for some z' G [z,z + 5.20]: 

°° 1 . x9 . 1 2m - 1 / 1 \ 1 , , 

r^Jv{t)dt> > 1 — -. (110 

^' -2m Iz' -\ 2i// 7(z + 5.20) ^ ' 

Next, consider the case v = m + (1/2). We assume m > 1 and z > 2m + 1. Using ([I], eqn. 
11.3.36), we can write 



m—1 

oo 



z 



•^^ k=l 

^ m—1 

- 2^^{'^^'^^'''^ + «^m+(l/2)(^)' + 2 Jk+imi^f)- 

k=l 

Then, using the lemma from the previous section, we get the following, for some z' G [z,z + 5.20]: 

/■°° 1 . ^9 , 1 2m / 1 \ 1 , , 

/ r^JJtfdt> r> 1 r- (111 

^' -2m + l7z'-\ 21// 7(z + 5.20) ^ ^ 

We now proceed to lower-bound the probability of observing a fine-scale element a. The following 
bound holds for any n > 2 and any r] < 1/e^. 

Pr[a <r]]> [ \f{k)\^dk 

= / ^p-"J„/2(27r/,/3)2.5op"-idp 



oo 



n / t-^J^,2{tfdt 

J27r/3/(Ar?) 



oo 



>n t-^Jn/2itfdt, 

Jn/(er]) 

where in the last step we used the fact that A > 27r/3e/n. Then, by equation (jlOOh . and using the 
fact that n/(er/) > 2e > 5.43, we get 

Pr[a<?7] > n(l--)——^ > n(l - = ^(l - -]. (112) 

^ eri ' erj 
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Pr[a<7?] > ^(l-r)- 



n 



(113) 



C.4 The variance of h orthogonal to 6 

First, we give a simple upper bound on integrals of the form 

for A; > 1, z > 2z^ and v > 1/2. This follows from equation (I102|) : 



/ t-''J^{tfdt < / t 

J z J z 



°° fe 2 2 



dt 



V3 



r^-^dt 



{-i/k)r 



(114) 



We now use this to upper-bound the integral 

/ Fo(r)V"-'^dr, 



'l/iX-qe) 

for /c > 1 and 7] < l/2e^. We write the following: (for the last step, recall that A < 2 • 27r/3e/n, which 
implies 2-^(3/ (Xije) > n/{2rje^)) 



Fo{rfr^-'dr = — / J„/2(27r/3r)V-'=dr 

l/(Ar?e) '-'o Jl/{\r,e) 



n 



< 



J2iTl3/{Xr)e) 

n 4 /27]e^\k 
So 7r\/3fc ^ 



Jn/2{trt-''dt-{2TTP) 

■ (2^/3)'=-i. 



k-1 



(115) 



In a similar way, we can upper-bound the integral 

/•OO 

/ F^(r)V"~'=dr, 

Jl/iXve) 



l/(A»?e) 

for A; > 1 and ?] < ^^(l/2e^). First, note that 

^o(^) = -ytr-"/V(„/2)+i(27r/3r)(27r/?). (116) 

(To see this, write Fo{r) = yj^(27r/3)''/2g(27r/3r), where g{x) = x-"/V„/2(x). Then F^(r) = 

y/^(27r/3)"/2g'(27r/3r)(27r/3), where g'{x) = -x-"/V(„/2)+i(x), see [I] eqn. 9.1.30.) Then we write 
the following: (for the last step, recall the fact that A < 2 • 27r/3e/n, which implies 27r/3/(A7/e) > 
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n/(2r/e2)) 



oo roo 



/{Xrie) 



So Jl/{X'rje) 



./(n/2)+i(2vr/3r)^r-'=dr.(27r/3) 



n 



< 



So J2iTl3/{Xrje) 



5*0 irVsk V n 



(117) 



We now combine this with the results of section 3, to show a bound on the variance of b perpen- 
dicular to 9, conditioned on a < 77. 

We can simplify equations ([30]) . ([37|) and ([55|) as follows: 



iBrlsih- 



da 







,n+l 



<\{n-2) 



4^" -J I Foirfr^'-^dr, 

'l/(Ar?e) 



r r r r da 

< lArlAlh^^ < 



7^2 1 /"OO 2-17 pco 

a«+i-2(n-l) AA/(A^e) ^ 2(n - 1) A A/(A^e) 



0< / /c./ci/2-^<(2n + 9 + ^)eA / 



Fo{rfr''-'^dr. 

l/(Ar,e) 

Plugging in our bounds for the integrals on the right hand side, we get: 

25 {2r]f 



iBrlsih 







n+1 



< 



r r r r da 5 

0< / lArlAlh^< 

Jo 



■ a^+i - n - 1 V5( 



5*0 n 
^(2r/)2(27r/3)2 + 



(27r/?)2, 

3600 {2r]f 



'0 



S'o n2 



(2vr/3)' 



0<f/c./c./2^<(l + A)^^(W)^- 



o 



n — 3/ S'o n 



Substituting into equation (f26]) . and using ()113l) . we get: 



E{lF{I -ee')b\a<r]) 

5.20/ 1 
< 11 + 



n — 1 



P'(M0{2r,f + 40(2r?)2 + ^(2^)' + 50(2^)=^ + ^(2r?)^) (ns) 



n — 3 



< (1 + 0(i))/3%(3536 + 260(2r?) + 0(i)) • 4. 
Using our assumption that 27] < 1/e^, we can rewrite this as 



Fibril - ee^)h\ a<r])< r//32(14300 + 0(i)). 



(119) 
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C.5 The variance of b parallel to 6 



We also get a bound on the variance of b parallel to 6, conditioned on a < ij. 
Substituting into equations ([76l) . ([82]) and ([93|) . we get: 



da 

KBrKBlK2——-r 



2'l 

< -(2r?)3(2vr/3)2, 



< fKArKAiK,^^ < ^(2r?)(2vr/3)^ + ^ 
da 700 (27?)3 



5300 {2'qf 



(2vr/3)= 



0< / KcrKciK2- 



,n+l 



/o « 
Substituting into equation (175]) . and using (|113p . we get: 



(2vr/?)^ + ^(2,)3(l + ^ 
bo n bo \ n — 6 



(2vr/3)^ 



• I a < ??) < ^ (l + (l2(2r?) + 100(2r?)3(f + i + I + 5 + 

< (1 + 0(i))/32 (63 + 3120(27y)2 + 0(1)) • 2. 
Using our assumption that 2r] < 1/e^, we can rewrite this as 



25 ■ 
n— 3 • 



E{{b.9)^\a<v)<p'{242 + 0{^)). 



(120) 



(121) 



D Spherical Shells 

D.l Spherical shell with square cross-section 

In this section, we give a heuristic analysis of the spherical shell with square cross-section. 

We can write / = h — g, where h is C times the indicator function of a ball of radius P + 6 around 
the origin, and g is C times the indicator function of a ball of radius /? around the origin. Then its 
Fourier transform is f = h — g, where h and g are calculated as in Section 4. We can write this as: 

f{k) = Fom, 

^oip) = ^(^ + <5)"/'^n/2(2vr(/3 + 5)p) - j„/2(2vr/3p). (122) 

We are interested in the case where 5 <C /3. Note that interference between the two Bessel 
functions begins to play a major role when p > 1/{2-it5). We claim that Fq{p) decays quite slowly, 
out to distance p ~ 1/{2it5). 

It will be convenient to define 

i^(/3)=/3-/2j„/2(27r/3p), (123) 

so we have 

Fo{p) = ^{K{P + 5)-Km. (124) 
We can approximate Fo{p) by a simpler expression. First, when 6 < (3/n, we can write C as follows: 

^ ' =M^. (125) 
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Also, when 6 is sufficiently small (we will elaborate on this point later), 

^O(P) « -^^K'iP). (126) 

A straightforward calculation (see [T], equation 9.1.30) shows that 

K'iP) = (27rp)/3-/2j(„/2)-i(2^/3p). (127) 

Note that K'((3) is roughly 2irp times larger than K((3), so we expect the approximation to be 
accurate when 6 < l/(27rp), or equivalently, when p < l/(27r5). Combining the above equations, we 
get the following approximation for Fq{p): 

^»<^> - ^i7m ■ -k-' ■ (^"')^"'^^<»/..-i(2-/'rt ^^^^^ 

Yo p{n/2)-i '^(n/2)-i(27r/3p) when p < l/(27r5), 

where Sq is the surface area of the sphere in M" (note that Bq = So/n). 

Compared to the Fourier transform of the ball (Section 4), this function decays more slowly, out 
to distance p ~ l/(27r5). Thus, when we apply the curvelet transform, with significant probability, 
we can observe fine-scale elements a < r], where r] shrinks proportional to 5. This suggests that a 
very thin spherical shell (i.e., 6 very small) allows us to find the center with very high precision, 
proportional to S. 

D.2 Spherical shell with Gaussian cross-section 

In the next few sections, we will prove Theorem IH for a spherical shell with Gaussian cross-section. 

We begin by proving upper and lower bounds on the normalization factor Cf. The following 
identity will be useful: (this follows from the definition of Fq(p) and a change of variables) 



r Fo{pfp^-'dp = Cje^p'^-^ I " " ^^^[-^eh'')J^^,,i^)^,{t)Hdt. (129) 
Then the norm of f{k) is given by: 

1^ poo 

/ \f{k)\^dk = So / Foipfp'^-'dp 

JR" Jo 

Jo 

= 5oC|e"/32"-2(Ar, + iV2), 
where we split the integral into two parts, 

Ni = / expi-^eH^)J^^/2yiit)hdt, (131) 
Jo 

POO 

/ exp{-^eH^)J^^/2)-iitftdt. (132) 

Jn-2 



N' 
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We now prove upper bounds on A'^i and For A'^i, using trivial upper bounds on exp(— x^) 
and Jy{xY (see [1], eqn. 9.1.60), we get: 



Ni < / ^tdt = i(n - 2)2. 
Jo 

For N2, using the upper bound on Jy{x)'^ from equation (|1U2|) . we get: 

/>oo 
Jn—2 



(133) 



4 -/2^1 



(134) 



Substituting into (|130p . we get 

\f{k)\^dk < 5oC|e"/32"-2(i (n - 2f + iJ^l). 



We assumed e < j^f^zfj^i ^^'^ it is easy to check that this implies j{n — 2)^ < 2y |i. Thus 

|/(^)|2dfc < 5oC2e"/32'^-2 . 



21 

3 £- 



Setting the left side equal to 1 implies a lower bound on C 



(^2 y ^ 1 fs 



Next we show lower bounds on A^i and For A^i we have a trivial lower bound 

Ni > 0. 



(135) 



(136) 



For N2, we use the lower bound J,/{x)^ > cos'^{9,y{x)) (see equations ()103p and (jlOip ). For 



convenience, we define 6{x) = 9^{x), suppressing the u subscript. We get: 

2 

N2>- exp{-^e^t'^) cos'^{eit))dt. 

Jn-2 



(137) 



Note that 9(t) is a monotone increasing function (equation (jl06p ). hence it is one-to-one and has 
a well-defined inverse. We make a change of variables, r = 0{t), t = 9~^{t): 



N2>- eM-^e\0~\T)f)cos^r 

^ Je{n-2) 



T)dT. 



Note that, whenever r = 9{t), we have {9 ^)'{t) = gr^- Hence, by equation p06p . 



1< (0-1)'(t) < ^, for r>0(n- 2). 



(138) 



(139) 
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Also note that 

e'\T) = e-\e{n - 2)) + f {9-^)'{x)dx 

Je{n-2) 

<n-2 + ^{T-e{n-2)). 

Substituting in, we get: 



(n-2) (140) 

2 



2 

N2>- exp(-^e2 . („ _ 2 + ^(r - e{n - 2))f) cos^ rdT. (141) 

7!" 7e(n-2) 



We will use the following simple fact: if a function / is nonnegative and monotone decreasing on 
the interval [a, oo), then 

/•oo roo 

/ f{x)cos^xdx>\ I f{x)dx. (142) 

This follows because 



/•oo °° ra+{k+l)7r 

/ f {x) cos^ xdx = / f{x)cos^xdx 

J a j(j_Q a+fcTT 

>^f{a + (k + l)7r) / cos^ xdx 

k=0 Ja+k7z 
oo 

= ^/(a + (A; + l)7r)f 
k=0 

oo /.a+(fc+2)7r 

= V/(a + (A; + l)7r)i / 

fc=0 Va+(fc+l)7r 
°o /.a+(fc+2)7r 

>^i/ f{x)dx 



(143) 



fe=0 

fOO 



/•oo 

= 5 / /(^)'^^ 

Using the above fact, and a change of variables, we get 

1 /•oo 

N2>- exp(-2^£2 . _ 2 + 2 _ ^(^ _ 2)))2)^^ 

Vr J6l(n-2)+7r 

^ r ^ 1 2 2x , 

= -— / exp(— X )dx 

V3 1 r / 2, , 

exp(-y )dy. 



(144) 



Recall that we assumed £ < This implies -^eln — 2 + ■^) < So, substituting in and 

"i"-^ v27r ^ VJ v27r 

integrating numerically, we get that 



iV2 > exp(-y2)dy > 1. (145) 

V27r£7 1 4e 
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Substituting into (|130p . we get that 

[ \f{k)fdk > 5oC|e"/32«-2(o + ^) = SoCP'^-'P 
Setting the left side equal to 1 implies an upper bound on C|: 



2^n-l,32n-2 1 
■ 4- 



(146) 



^ - 5oe"-i/32"-2- 



(147) 



D.3 The low- frequency components 



Next, we show that f{k) has very little power at low frequencies, corresponding to coarse scales 
o > 1. This justifies our use of the curvelet transform, which effectively ignores these low-frequency 
components (recall Theorems [T] and [2]) . 

The total amount of power at frequencies less than z (for any z > 0) is given by: 



\f{k)\^dk = So / Foipfp'^-'dp 

\k\<z Jo 

2TTf3z 



Jo 



(148) 



(We used equation (jl29p .) Using a trivial upper bound exp(— x^) < 1, and the upper bound for 
Ju{x) (when x is small) from equation (|96]) . we get: 



•<2^n p2n~2 



2ttI3z 



Ifik^dk < SoCje'^P 

\k\<z JO 

C /^2 n2n—2 
= ^O^fS P 



■Kin - 2) Vn - 2 



te \"-2 



tdt 



^n-2 



(149) 



7r(n - 2)"-i n 
Using our upper bound on C| (equation (jl47p ). we get 



(27r/?z) 



\f{k)?dk < 4e • -— ^ • -{2^Pzr 

\<z 7r(n-2)" i n 



< 



Ae 



(150) 



7re2 (n - 2)' 



• (27r/3z)". 



Now, fix z 



n-2 
27r/3e- 



Recall that A > hence l/\<z. Then we have 



\f{k)\''dk< 



|fe|<l/A 



Ae e 
-^5- 



vre 



(151) 



Recall that we assumed e < So the frequencies below 1/A only constitute a small fraction of 

the total probability mass. This justifies our use of the curvelet transform, with this choice of the 
parameter A. 
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D.4 The probability of measuring a fine-scale element 

We give a lower bound on the probability of measuring the scale variable to be small, a < rjc, where 

6 (n-2) , , 

We will show that a < ijc with at least constant probability. 

First, we write Pr[a < r]c] as follows, using Section 3.1 and equation (jl29p : 



Pr[a < r],] > Sq / Fo{pfp'''^dp 

Jl/(\Vc) 

/•oo 

J27T0/(XVr) 



(153) 



'27r/3/(A»7c) 

The lower limit of integration can be simplified, by substituting in the definitions of A and ijc, 

- 2.^(M£)-V£(!i^)-' = 2.«2,i)-' = f = i (154) 



Xr)c \n — 2/\^e/ 6 e 

This integral is similar to the integral which we encountered earlier (equation ()132p ). We can get a 
lower bound using the same technique (equations (jl37p - ()144p l: in particular, note that 1/e > n — 2, 
as required in that calculation (this holds because we assumed e < l/(n + 2)). This leads to: 

Pr[a < 77c] > 5oC|e"/32"-2 . T eM-y')dy. (155) 



The lower limit of integration can be written as 



1 /I 27r\ 1 / 27r \ 1 2 

e - + ^] =^ l + ^e] <^-2 = x (156) 



V2^ Ve V2^\ VS / V2^ V vr 

where the last inequality follows because e < i. Thus we can lower-bound our integral as follows: 

Pr[a < r?J > SoCje^P^''-' ■ / exp(-y2)dy 



2.n«2n-2 0.15 (157) 



e 

y2^n-l n2n-2 



= (0.15)5oC|e"-V" 
Now, using our lower bound for C? (equation (jl35p ). we get: 



Pr[a < Tjc] > (0.15) • W| > 0.045. 



(158) 



D.5 Some more integrals 



Our next goal is to bound the variance of b. We begin by proving upper bounds on certain integrals 
involving -Fo(^) and i^o(r). Then, in the following two sections, we will bound the variance of b in 
the directions orthogonal and parallel to 9. 
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First, we consider integrals of the following form, where k > 1: 



oo 

Fo(r)V"-'=(ir. (159) 



Using the definition of Fq^t), and a change of variables 

\2n-kr^_^2^n/r,\k-ln2n+k-3 / 1^2.2 



dt 



We will upper-bound this integral, assuming that a > First, using equation ()102p . we have 
J{n/2)-iitf < and we get: 

Next, we use a simple inequality: t^^^ > {211^0)^^^ , whenever t > 27r(3a. Thus, 

J a TrVSa" ^ J2npa 

The integral on the right hand side can be bounded as follows: 

/ eKp{-^e'^f)dt = / exp(-r^)(ir • V2^- < \- ■ V2^- = 

Substituting in, we get: 



/ Fo(r)V"-'=dr < C7|£"-i/32"-2 . —^3^, for a > 

J a V D ^ 



n-2 
27r/3 ■ 



Next, we consider integrals of the following form, where k = 1,2: 



(160) 



00 

F^(r)V"-^dr. (161) 

In order to calculate -Fo(r), we write -Fo(^) in the following form: 

Fo(r) = C^e"/2(2vr)/3" • (27r/3)("/2)-i . ^(r) . K(27r/3r), 

where ^ 

i?(j;) = exp(-7r£^/3^a;^), K(x) = ^(„/2)-i '^(n/2)-i(a;)- 

Then FQ{r) can be written as 

F^(r)=F('^)(r)+F(')(r), 

where 

F(S'^)(r) = Cye"/2(27r)/?" • (27r/3)("/2)-i . _?^'(^) . K{2Ti(5r), 
F^^\r) = Cy£"/2(27r)/3" • (27r/3)("/2)-i . . K' {2^(ir){2^p). 
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\j.{nl2)-l-^ri/2{x) 



We can expand this out. Note that (see [T], equation (9.1.30)) 

H'{x) = eicp{-TTe^l3'^x'^){-2TTe^l3'^x), K'{x) -- 

Substituting in, we get 

F^''\r) = C7^e"/2(2vr)/?"-^^^ • exp(-7r£2/52r2)(-2^e2/3V) J(„/2)_i(27r/3r), 
Fil'\r) = C^e«/2(2^)/3"-^-l-^ . exp(-7r£2/3V2)(-l) J„/2(2vr/3r)(27r/?). 

Thus -Fo(r) is given by: 

Fi{r) = Cfe^'\2T,)(3^-^^^eM-'^e^P\^)\-27,e^P^ (162) 
Substituting into our integral, and performing a change of variables, we get: 
/ F^(r)V-'^dr = C|e"(27r)^-i/32"+'=-3 / exp(-^e2^2) (^^2^^^^^^_^(^) ^ 27r/3J„/2(t) ) 

Ja J2TTl3a ^ ^ 

(163) 

We will upper-bound this integral, assuming that a > First, using equation p02p . we have 
Jn/2{^f' < ^"^^ similarly for J(„/2)-i(i)^- So we get: 



F^(r)V"-^dr < C72e"(27r)'=-V^"+^"^ 



exp(-^e2f2)/^2^^_^27r/3 



Changing variables and rearranging, we get: 
F^(r)V"-'=dr 

2n(^\k~l a2n+k~3 



2iT£l3a 



exp(-r^)(V27re/3r + 27r/?) 



2^J^ ^ fc_2 



r'2^n(r, \k-l o2n+k-3 / l_c-'i^-2 



exp(-r2)(27re2/32^2 ^ 2(27r)3/2e/?2^ + (2vr)2/32)^ 



C^e"+^-2V2^V'"+'-' • ^ / _ exp(-r2)(2^eV2 ^ 2(2^)3/2er + (27r)^ 



27re/3a 



rk-l • 



We can handle integrals of the form 



roo 

/ exp(- 



dT 



1 ^ 

■ 2 • 



as follows. When ^ > 0, we have 

[ exp(-r2)^ ^ -J f exp(-r2)(ir < 

— 1, we have 

/■OO f'OO 

/ exp(— r^)r(ir < / exp(— T^)r(iT = — ^ exp(— r^) 

Ja Jo 



When 



(165) 



(166) 



(167) 



(168) 
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When 



-2, we have 



\ exp(-r2)r2dr < / exp(-r2)T2(ir = -i exp(-r2)r -/ -i exp(-r2)dr = ^ < 0.45. 

(169) 

Now, we can upper-bound our integral in the k = 2 case: 



dr 



exp(— r )(27reV2 + 2(27r)3/2er + (27r)2) 



C|e"(2.)/32"+^ . -I (ie^ + V2vr. + -^^) , for a > 



^2 £/3c 



27re/3a 2 ^ 

n 

2^' 



We can also upper-bound our integral in the k = 1 case: 



/oo 
exp(-r2)(27reV2 + 2(27r)3/2er + {2Trf )dT 

< ^^£"-1^2^/32" . (^2^e2 . ^ ^ 2(27r)3/2e . i + (27r)2 • ^ 



^^£"-1^2^/32" • ^ (^e^ + \/2^e + vr^/^) , for a > 



n 



27r/3 



(170) 



(171) 



D.6 The variance of b orthogonal to 6 

We now bound the variance of b orthogonal to 6, conditioned on observing a < ijc- Recall that 



Vc 



_ S (n-2) 

We start with the results of Section 3. From equation (j30p . we get: 



poo 

<\{n-2)\rie j FQ{rfr''-^dr. 
h/{\ne) 



(We used the fact that 1 < r • Ar/e, for all r in this interval.) From equations (j37p . 
get: 







< / lArlAih^, < • y / F^(r)V"-2dr+ 

2(n-l) XJi/l^Xrie) 
^2 />oo 

24e2Ar/2 / FQ{rfr''-^dr. 

Jl/iXve) 



(172) 
and we 

(173) 



2(n - 1) 

(Again, we used the fact that 1 < r • Xrje, for all r in this interval.) From equation (|55|) . we get: 

/■^ dn 

0< / /cr/ci/2— T<2(n + 10)eA / Fo(r)V-2dr. (174) 
(We used the fact that | 2(n - 2)(1 + ^) = | 2(n + 3 -h ^) < 2(n 10), assuming n > 4.) 
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Now we fix = ?7c, and we upper-bound these integrals, using equations (|160p and (jl70p . (Note 



5 in-2) 



that, in order to apply these results, we must have > Recall that A = and r/c 

Also, recall that we assumed e < —. Then Xrin = 27r5 = 27re/3, and = o \ > tt^, as desired.) 
After some tedious calculation, we get: 



< (n - 2) • C|e"- V' 



^2 n-1^2n-2 _o2/i r,nnn_2 



do 

< / lArlAih—T < 5 • • e/3"(12000e" + 9e + 80), 

Jo 



< Icrlcih-^ < C}e--'p'-~' ■ 2600e/?2(l + ;i^) • 5. 

Then, by combining these equations and using our assumption that e < O(^), we get: 
rvc , \ da 

< C|e"- ^^"-2 . £^2^600006^ + Abe + 400 + 2{n - 2)120e + 2600(1 + :^)S) 



< C72e"-i/32'^-2 . e^2^4QQ ^ 2600 + 0(i)) • S. 



Next, recall from equation (jl57p that: 



Pr[a < > (0.15)5oC?e"-V2n-2^ 



(175) 
(176) 
(177) 



(178) 



(179) 



Finally, by substituting into equation (j26p and simplifying, we get a bound on the variance of 6, in 
the subspace orthogonal to 9, conditioned on observing a < ijc- 



E(?^(I 



a<Vc)< (n - 1)6/32(507 + 0(i))- 5. 



(180) 



D.7 The variance of b parallel to 9 

We now bound the variance of b parallel to 6, conditioned on observing a < r]c- Recall that r]c 

S (n-2) 

We start with the results of Section 3. From equation (j76l) . we get: 



da 

KBrKBlK2—^ 



< \{n - l)(n - 2) / Fo(r)V"-3dr. 



From equations ([82]) . (p3|) and ([85|) . we get: 

< / KArKAiK2—r < 2 / F^(r)V"-Mr+ 



l/(Ar,e) 

/>oo 

2 • fvr^ / Fo(r)V"-3(ir. 

Jl/CArje) 



From equation ([93|). we get: 



< / KcrKciK2—^ < 4(n2 + 3n + 3) / Fo(r)V-3cir. 



(181) 



(182) 



(183) 
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(We used the fact that ^(n - 1) + ^(n - 2)^(1 + ^) = zj(n2 + |n - | + < ^(n^ + Sn + S), 
assuming n > 4.) 

Now we fix = rjc, and we upper-bound these integrals, using equations (jl60p and (jl7ip . (Note 

that, in order to apply these results, we must have > Recall that A = and r/c = ~ • 

Also, recall that we assumed e < Then Xrjc = 2t:5 = 27re/3, and = ^^^^^ > as desired.) 
After some tedious calculation, we get: 



^'^ da 

KBrKBlK2 



n+1 

da 



< C2e"-i/32«~2.n2.i20e2/32, 



< / KArKAiK2—T < C7^e""^/3^"-^ • /3^(25200e^ + 60e + 132), 



da 



< / KcrKciK2—T < C^e""^/?^"-^ • 2400(n^ + 3n + 3)e^p\ 
Jo 0- 

Then, by combining these equations and using our assumption that e < O(^), we get: 



(184) 
(185) 
(186) 



Vc 



KArKAlK2 - 2KBrKBlK2 + KcrKciK2 



da 

-,n+l 



(187) 



< C|e""^/32""2 • /32(25200e^ + 60e + 132 + 2 • 12idr?e^ + 2400(n2 + 3n + 3)6^) 

< C72e""i/32n-2 . ^2(^32 ^ 0(;^)). 

Next, recall from equation (jl57p that: 

Pr[a < ?7c] > (0.15)5oC|e"-i/32"-2. (188) 

Finally, by substituting into equation (I75p and simplifying, we get a bound on the variance of 6, in 
the direction 0, conditioned on observing a < rjc- 



E{{b-9f\a<nc)<P\23 + Oi^)). 



(189) 



E A Fast Quantum Curvelet Transform 
E.l The Discrete Curvelet Transform 

First, we argue that the discrete Fourier transform approximates the continuous Fourier transform, 
in the sense described in Section 6.1. This follows from the definitions of the different transforms. 
Recall the continuous Fourier transform that takes a function on to a function on R": 

:Fcont{fm = I /(x)e-2-*^-^^dx, (190) 
^FcoUam = [ g{k)e'-''-'dk. (191) 

JR" 

Now consider the Fourier transform that takes a function on the cube C = [—L, L)" (or equivalently, 
a function on that is periodic with respect to the lattice (2LZ)"') to a function on the (dual) 
lattice C = (g^Z)". We refer to this as the "semi-discrete" Fourier transform: 

:Fsen.^{f){k) = {^T'^ [ f {x)e~^-^'^-^ dx, (192) 

Jc 

^^eUam = E 9{k)e'''''-'- (193) 
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Also recall the discrete Fourier transform, that takes a function on Z = (crZ)" n [—L, L)" to a 
function onZ = {^Zrn[-l,^y-: 

:FUfm = (194) 
^Ms{9){S) = i^r/' g{k)e'-''-'. (195) 

We are given a function fcont on M" that vanishes outside the cube C. We define a function fsemi 
on C by restriction, fsemi = fcont\c- Then it follows from the definitions that fsemi = {^T^"^ fcont\c- 

Recall that fcont has most of its probability mass inside the cube 2^)"- Then the same 

should be true for fsemi- Now define a function fc^is on Z by restriction, f^is = fsemilz- Then, using 
the definitions, we see that fdis ~ f^^^/semiU- 



Note that an example of a discrete curvelet transform (over or M^) can be found in [6]. There, 
the frequency space is partitioned into concentric cubes according to the scale a, and these are 
divided into wedges according to the direction 6. For our purposes, we will use a tiling based on 
concentric balls (in M"), which more closely approximates the continuous curvelet transform defined 
in Section 2. 

As a classical computation, the discrete curvelet transform can be implemented using the fast 
Fourier transform [6j (see in particular the "wrapping" method). We will use these ideas to implement 
a quantum curvelet transform. The following discussion will be self-contained; but for readers who 
are familiar with [6], we mention that we omit the "wrapping" step. Our transform produces curvelet 
coefficients that are somewhat oversampled, but this does not cause any problems in our situation. 

E.2 Constructing the Window Functions Xae(^) 

We will construct two families of window functions Xagi^)^ which the operation X can be 
performed efficiently. 

First, suppose we have some partition of the frequency domain into disjoint subsets, (Zm)" = 
Uae'^aS' such that given any point k £ (Zm)"", we can efficiently compute which set S^^ contains 
k. Define the window functions to be the indicator functions for these sets, 

I 0, otherwise. 

Then the operation X can be implemented efficiently: it simply maps |/c)|0,0) \k)\a,9), where a 
and 9 denote the set S s that contains k. 

Unfortunately, these window functions are sharply discontinuous, so the resulting curvelets are 
not very well-localized in space. This makes them poorly suited for the applications proposed in this 
paper (recall that the results of Sections 3, 4 and 5 required window functions that were C^-smooth). 

Smooth window functions are more challenging to implement, because the supports of the func- 
tions Xa gi^) necessarily overlap. Thus, at a given point k, the operation X must create a super- 
position of many values of a and 9. These superpositions can be complicated: for instance, if we 
imagine that the tiling of frequency space looks (locally) like an array of n-dimensional cubes, then 
a significant amount of volume lies near the corners of the cubes, and each corner point touches 2" 
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different cubes, so we would have to prepare superpositions of 2" different values of a and 6. This 
seems impossible for many choices of the window functions. 

However, the above example also suggests a solution to the problem. We can use spherical 
coordinates, which look locally like Cartesian coordinates (except at the poles). If we define the 
window functions to be products of simpler functions, each depending on a single variable, then we 
can prepare these superpositions efficiently. We now demonstrate this construction. 

First, recall the definition of spherical coordinates in M": we have {r,(pi, . . . , (pn^i), where r E 
[0,00), (pi, . . . ,(f)n-2 e [0,7r]U{"undef"}, and € (-vr, vr] U {"undef"}. We use the value "undef" 



to represent points on the "poles" of the sphere, e.g., if (pj = or tt, then (pj+i = ■ • • = (pn-i = 
"undef" (a similar situation arises when r = 0). 

Cartesian coordinates are written in terms of spherical coordinates as follows: 

xi = rcos(f)i (or if undefined), (197) 

Xj = r sin ^1 • • ■ sin (f)j-i cos cpj (or if undefined) {j = 2, . . . ,n — 1), (198) 

Xji = r sin (pi ■ ■ ■ sin (f>n—i (or if undefined). (199) 

The reverse mapping is given by: 

r = ^xl + ---+xl, (200) 

01 = arccos(xi/r) (or "undef" if r = 0), (201) 



(f)j = arccos(xj/ (r sin (pi - ■ ■ sin (pj-i)) (or "undef" if (pj-i G {0, tt, "undef" }) (j = 2, . . . , n — 2) , 

(202) 

(pn-i = sign(a;„) arccos(a;„_i/(r sin(^i • • •sin(^„_2)) (or "undef" if (pn-2 € {0,7r, "undef"}). (203) 

Next we will define discrete values for the scale variable a and the direction variable 9. In the 
notation, it will be convenient to represent the scale variable by s instead of a, where 

a = 2-^ (204) 

We will then define window functions Xg gik)- These will be products of radial and angular compo- 
nents (we write k = (r, (p) using spherical coordinates) : 

X,/k)=WsiXrX/$). (205) 

Here, A is a parameter that sets the radial scaling. For future use, we define the function c : [0, 00) 
M, 

, , Icosx, < X < 7r/2, 
1^0, X > 7r/2. 

We begin with the scale variable a = 2~^. We fix the cutoff values Smin,Smax € Z, where 
1 < Smin < Smax- Then we let s G {smin, Smin + 1, ■ ■ ■ , Smax} U {"coarse" , "fine"}. 
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Figure 2: A family of radial window functions: w^c 

We define radial window functions Ws{r) as follows: 

'c(f (2^ - r)/2*-i), 2^-1 <r< 2 



',wi,W2 and tyw- 



Ws{r) 



,(r) = < 



S 

1 

c(f(r- 2*)/2^), 2*<r<2*+\ 

0, otherwise, 

1, < r < 2'*'"^""^ 
^0, r > 2*™'", 



0, 



< r < 2'*™'^==, 



An example is shown in Figure [2j It is easy to check that 

"^Wsirf = l (Vr>0). 



(207) 



(208) 



(209) 



(210) 



(Note that at any given point r, at most two of the functions Ws{r) are nonzero.) 

We let the direction variable 9 take on values in the set Gs(S'"^^). (Assume for the time being 
that s ^ {"coarse", "fine"}; we will handle those special cases later.) The set Gs{S^~^) contains 
grid points on the sphere 5""^, defined using spherical coordinates, with angular spacing 7r/2l^*/^^ ~ 
7r/\/2^ = vry^. This set is defined recursively: 

G,{S^) = {^t/2r«/2l \ teZ, 0<t<2- 2r'^/2l _ (211) 
GsiS'') = {7rt/2r^/2l I i g 1 < t < 2^^/21 _ i| ^ GsiS''-^) U {0,^} x {"undef"} {k > 2). (212) 

(See Figure [3] for an example.) 

We define angular window functions v Jcj)) as follows: 



n-l 



where 



= n '^s,e,{(t>j), 
Us,e,{cP,) = c{2^'^/^^\<l>,-ej\/2). 



(213) 



(214) 
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Figure 3: Gq(S'^), the set of grid points on the sphere S*^ C M^, with angular spacing vr/S. 



This requires some further explanation. 

Intuitively, Ugfiji^j) is a one-dimensional "bump" function centered around 9j, and v^^cj)) is a 
product of these functions. 

Note that, for the first n — 2 coordinates j = 1, ... ,n — 2, Us^Oji'Pj) is defined on the interval 
[0,7r], whereas for the last coordinate j = n — 1, Usfi^{(t>j) is defined on the circle (— tt, tt]; in this 
latter case, we interpret — Oj\ as the shortest-path distance around the circle. 

Also, in the definition of v^^{(f)), we simply omit those factors that have dj = "undef" or <pj = 
"undcf". We claim that this is a natural thing to do, in that it yields a simple geometric picture. 
Intuitively, Oj = "undef" means that 6 is located on a pole of the sphere, with some coordinate 9i 
{i < j) equal to or tt. Then this construction produces a bump function that covers a circular 
region around the pole, and so does not depend on (pj. On the other hand, if 4>j = "undef", then (j) 
is located on a pole of the sphere, with some coordinate (pi {i < j) equal to or tt. If 9i ^ then 
'Usfii{4'i) = 0, hence v^^{(p) = 0, independent of (pj. If 9i = (pi, then 6 is located on a pole, hence 

ff{(p) does not depend on (pj . 

Note that at least the first {j = 1) factor will always be defined, since 9i is always defined 
whenever 9 G Gs{S^~^), and (pi is always defined whenever k (we can ignore the case of A: = 0, 
because it is relevant only when s = "coarse" , in which case we will not use these angular windows) . 

Finally, in the special case where s = "coarse" or "fine", we do not resolve any directions 9. 
Instead, we fix = "undef", and we define the angular window to be trivial, ^<p) = 1. 

We now show how to perform the operation X that maps 

|fe)|0,0)^|^)^X,/(fc)l^,^"). (215) 

sfi 

We will do this by converting k to spherical coordinates {r,(p), performing an operation X' that 
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creates the superposition over s and 9, then converting back to Cartesian coordinates: 

|^)|0,0)|0,0) ^ |fc)|r,(^)|0,0) 

^\k)\r,$)Y,_^Ws{Xr)v^g{$)\s,e) 

s,e 

= \k)\rJ)^^Xs,e<k)\s,e) (216) 

s,9 

^\k)\O,6)^^Xs,e<k)\s,0). 

sfi 

The operation X' is implemented recursively, acting on the variables r, ^i, . . . , </!>n-i one at a time: 
/ / The "s" register 

If Ar < 2*™*", then let si = "coarse" and S2 = Smin- 

Else if Ar > 2*'"°'^, then let si = Smax and S2 = "fine". 

Else, let si = [lg(Ar)J and S2 = [Ig(Ar)]. 

If ,si = S2, then set the ",s" register to 

Else, set the "s" register to Ws^{Xr)\si) + u^^j (Ar)|s2). 

// The "Oi" register 

If .s € {"coarse", "fine"}, then set the "^i" register to |"undef"). 
Else, begin: 

Let n = (7r/2-)L0i(2V7r)J and T2 = (7r/2-) [0i(2V7r)l . 

If Ti = T2, then set the "^i" register to |ti). 

Else, set the "^i" register to ^is,ri + ^s,r2 I''"2)- 

/ / Note, 7^ "undef" , since otherwise we would have s = "coarse" 

// Note, if (pi G {0, tt}, then ri = T2, hence 9i G {0, tt} 

End. 

Recurse on the "^2" register. 

// The "0/' register, for £ = 2, . . . , n - 2 

If 9i_i G {0,7r, "undef"}, then set the "6'/' register to | "undef"). 

Else, begin: 

Let n = (7r/20L'^^(277r)J and T2 = (7r/2^) [0^(2771)] . 

If Ti = r2, then set the "^£" register to |ri). 

Else, set the "6*^" register to 'Us,^i((^£)|ti) + ^5,^2(^^)1x2). 

/ / Note, 7^ "undef" , since otherwise we would have, in some previous iteration A;, 
// (pk e {0,7r}, hence 9k G {0,7r}, and 9e-i G {0,7r, "undef"} 
/ / Note, if (p£ G {0, tt}, then ti = T2, hence 9^ G {0, tt} 

End. 

Recurse on the "0^+i" register. 
// The "6l„_i" register 

If 9n-2 G {0,7r, "undef"}, then set the "6'„_i" register to | "undef"). 

Else, begin: 

Let n = (7r/27L<^„_i(277r)J and T2 = {7r/2')\4>n-i{2' M]. 
If Ti = T2, then set the "^„-i" register to |ri). 
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Else, set the "9n-i" register to Us,Ti{4>n-i)\n) + tts.ra (</'n-i)|T2). 

// Note, (pn--! 7^ "undef", since otherwise we would have, in some previous iteration k, 
II (f)k G {0,7r}, hence 6k G {0,7r}, and 6n-2 G {0,7r, "undef"} 

End. 



This construction yields a fast quantum curvelet transform using smooth window functions. 
Note that we can carry out this construction using other choices of the function c(x), which lead to 
different window functions Xs e^^)- We only need c{x) to satisfy the identity c(x)^ + — x)^ = 1 
(for all < X < 7r/2). 

For instance, we can define 

/cos(/i(x)), 0<x<7r/2, 
4x) = < (217 
[0, X > 7r/2, 

where h{x) is any increasing function that satisfies /i(0) = 0, /i(f ) = f , and /i(-| — x) = | — h{x) 
(for ah < X < 7r/2). 

In particular, if we set h{x) = ^ sin^ x, then the resulting function c(x) is C^-smooth. Thus we 
get window functions Xg g{k) that are C^-smooth, and are qualitatively similar to the ones used in 
Sections 3-5 of this paper. 



F Quantum Algorithms using the Curvelet Transform 

F.l Single-shot measurement of a quantum-sample state 

Here we analyze a "continuous" analogue of our single-shot measurement procedure for finding the 
center of a ball. 

First, we claim that a < rj with probability > By Theorem [3l we have 

Pr[a < r?] > (0.19)r?(l - i). (218) 

Note that r/ > 0(i/^), which follows from the definition of t] and the fact that = uf] > vl3l2. This 
shows the claim. 

From this point on, all probabilities are conditioned on having a < ij. 

Without loss of generality, assume c = 0. The algorithm succeeds when it outputs a point close 
to 0. Let Hi be the projector onto the subspace orthogonal to 6, and let 112 be the projector onto 
the direction 9. We will show that, with constant probability, |ni6| is small and |n26| is not too 
large. 

Let X = |ni5|2, nx = E(X), and let Y = |n26p, fiy = E(y). By Markov's inequality, 

Pr[X > 3fix] < I, Pr[y > Sfiy] < i 

Then, the union bound implies: 

Fi[X < 3fix and Y < Sfiy] > 1 - Pr[X > Sfix] - Pr[Y > Sfiy] > |. 

So, with probability > 1/3, we have X < 3 fix and Y < Sfiy- 

We now rewrite this in terms of |ni6| and |n26|. From Theorem [3l we know that 

fix < r?/3'(14300 + (Qi/n)), fiy < (3\2A2 + (^2/^)), 
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for some constants Qi and Q2- So, we have 



|ni6| < \/3V^/3\/l4300 + (Qi/n), |n26| < ^3^^/242 + (Qz/n). (219) 

This shows that |ni6| (the error orthogonal to 9) is small. Indeed, substituting in our choice of 
r/, and using the fact that (3/(3 < 1, we see that 

|ni6| < ^/x. (220) 

However, |n26| (the error parallel to 9) is not so small. So the algorithm tries to guess this error 
and output a corrected point. It succeeds when 

\li2b- n(^/3V242 + {Q2/n))9\ < -j^fi. (221) 

Call this event E. The probability of E is the probability that a random point chosen uniformly 
from the interval ^/3(^^/242^\^XQ2/n} ■ [—1, 1] lies within distance -^/x of some fixed point in the 

(possibly smaller) interval \/3/3\/242 + {Q2/n) ■ [—1, 1]. This probability is lower-bounded by 

1 

Pv\E] > , ^ = ^ ^, (222) 

2^/3/3^242 + {Q2/n) 2^6^242 + (Qs/n) /3 

which is > f^(i^), since ^ = 1^(3 > v(3/2. 

So the algorithm succeeds with probability > $7(z^^). This proves Theorem [5l □ 

Finally, we argue that, when the grid G = (crZ)" Pi [—L, L)" is chosen properly, the discrete 
algorithm will behave like the continuous one. Recall from Section 6 that we can approximate a 
function / on M" with a function /2 on G, provided that most of the probability mass of / lies within 
distance L of the origin, and most of the probability mass of / lies within distance ^ of the origin. 
This holds for our algorithm, provided that: 

(The first condition follows immediately, since / is supported on a ball. The second condition follows 
from the decay of /; details omitted.) These conditions hold whenever 

L> R + i3, a< —€ — . (224) 

- - 600 /? 14300n + Qi ^ ^ 

Also, we argue that the discretization of the "direction" variable 9 will not introduce too much 
error in the output of the algorithm. Our algorithm constructs a line i = {b + X9 \ AG M}, and if 
9 were a continuous variable, this line would pass within distance 0{^/a(3) of the center c. Recall 
from Section 6 that the discrete curvelet transform resolves 9 within error ity^ in angular distance. 
Since h lies at distance 0{(3) from the center c, the error in 9 can increase the distance from ^ to c 
by at most 0{^/a(3). Thus the error in the output of the algorithm increases by at most a constant 
factor. 

We can bound the running time of our algorithm as follows. Let M be the number of grid points 
along one direction. Then M = 2L/a, and the running time is < poly(n, log M). Say we choose L 
and a so that the above inequalities are tight (up to constant factors). Then M < 0{R(3n/ fj,'^), and 
the running time is < poly(n, log R, log (3, log ^). 
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F.2 Quantum algorithm for finding the center of a radial function 

Here we analyze the "continuous" analogue of our algorithm for finding the center of a radial function. 
First, consider what happens for each i £ {1, 2}. 

We quantum-sample over a ball of radius R' around 0, then measure the value of /, and get a 
superposition over a shell of radius around c. If the shell is too large, it will not lie completely 
within our original ball, so we only get a fragment of the shell. However, if < R' — R, then we 
are guaranteed to get a complete shell. 

We claim that we observe /J*-*) < R! — R with constant probability. To see this, write: 

^ r ^1 ^1 volume of ball of radius R' — R around c 

Vi[(3^^> < R - R] = — 

volume of ball of radius R' around 

(Tit _ r>\n [ZZii) 

= lli— i^L = (1 _ i)" > e~i-2 > 0.30, 

using the fact that 1 — x > e~^^''^^^ for all < x < 1/4 (recall that we assumed n > 4). 

Our algorithm does not know the shell's true radius /J^*-* , so it uses as an estimate. 

We claim that < f3^^^ < (3/2)/?(*), with constant probability. Observe that 



Pr < (2/3) (i?' - R)] = '^Q^^™^ °f ^^^^ ^^^^^^ (2/3) (i?' - R) around c 

volume of ball of radius R' around 
((2/3)(i?'-i?))" 



(226) 



(2/3)"(l-i) 



So 



Pr[(2/3)(i2' -R)< <R' -R] = {1- (§)")(! - ^)" 

> (1 -0.20)(0.30) = 0.24, 



(227) 



using the fact that n > 4. 

Next, we claim that we observe a^*^ < r]^^\ with constant probability. This follows from Theorem 

m 

Pr[aW < > 0.045. (228) 
From this point on, we take probabilities conditioned on a^*^ < r]^^\ 

Without loss of generality, assume c = 0. Let H^*^ be the projector onto the subspace orthogonal 
to 6^^\ and let be the projector onto the direction 6^^\ 

We claim that |n^*^5(*)| is small, and |n2''fe'-*^| is of order /3, with constant probability. We use 
the same argument as in the previous section, together with Theorem [H We define e^*^ = 5/P^'^\ 
We get that, with probability > 1/3, 

|nP6«| < Vs^in - l)eW/3« ^761 + (Qi/n), |n^*)6«| < V3/3« V23 + (^2/^^), (229) 

for some constants Qi and Q2- 

Using the definition of e^^\ and the fact that Z?^*-* < R' — R = {n — l)i?, this implies that 

inP^W] < y3(n - 1)V5R^/W^T{Q^, \Jlfb'^^ \ < V2,{n - l)Ry^23T{Q^/^. (230) 

The algorithm carries out this procedure twice, for i = 1 and 2. With constant probability, this 
produces two lines, ii = {6^^) + A^(^) | A G M} and £2 = {6^^^ + Ai9(2) | A G R}, which both pass near 
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Figure 4: 



the point c. The algorithm then checks that these Unes are nearly orthogonal, and if they are, it 
returns the point on ii closest to ^2- A straightforward calculation shows that this point is given by 
^^(1) +6(1). 

First, we claim that the lines ii and £2 are nearly orthogonal (|^^^^ • ^^^^| < 3/4) with at least 
constant probability. 

We want to upper-bound the probability that \9^^^ ■6^'^^\ > 3/4. Recall that these are independent 
random vectors, chosen uniformly from the unit sphere S""^^ in M". It follows that 

■ > 3/4] = Pr[|xi| > 3/4], (231) 

where x = (xi, . . . random vector chosen uniformly from ^. Note that E(x) = 0, hence 

^{xi) = 0; also, E(|xp) = 1, hence E(a;f) = 1/n. Then by Markov's inequality, 

Pr[|xi| > 3/4] = Pr[x? > 9/16] < f ^ < | (for n > 4). (232) 

Thus, we observe • < 3/4^ with probability > 5/9. (This is a rather weak bound, especially 
when n is large, but it is adequate for our purposes. Actually, it is the case that |^(i)-^(2)| < 0(1/V^), 
with probability > f^(l).) 

Next, we claim that when ii and £2 are nearly orthogonal (l^^^-* • 9^'^^\ < 3/4), the point on ii 
closest to £2 (call it pi) is close to c. 

Let qi be the point on ii closest to c, and let pi be the point on £1 closest to £2- Similarly, let q2 
be the point on £2 closest to c, and let p2 be the point on £2 closest to £1. (See Figure [H) 

We know that qi and q2 are both close to c: |gi — c| < A, and 1(72 — 51 < A, where A = 
\/3(n — l)V6Ry/761 + (Qi/n). Furthermore, pi and p2 are close together: \pi — P2\ < 2A. 

Now suppose that pi is far from c: 

\pi-c\> 8A. (233) 

From this we will derive a contradiction. 
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First, note that pi is far from qi: 

\pi-qi\>7A. (234) 

Consider the hne £'2 = ^2 + (pi — ^2)- Define another point q'2 = Q2 + {pi — P2)- The hne £'2 is 
parallel to £2, it intersects £1 at pi, and it passes through q'2. Note that q2 is close to c: 1^2 — 
Note that pi is far from 

\pi - gal > 5A. (235) 
Note that qi and ^2 are close to each other: 

|gi - q'2\ < 4A. (236) 

We will use equations (I234p . ()235p and (I236p to show that the angle between £1 and £'2 is small. 
(See Figure H) We have o < 4A, 6 > 7A and c > 5A. 

Using the law of cosines, and the fact that x + ^ > 2 for all x > 0, we get: 

„2 I l2 „2 

+b -a 



26c 

1 / c 6 \ 
2\b^c ~ be) 
a" 



(237) 



> 1 

26c 



(4A) 



2 



3 

35 ^ 4- 



> 1 i '- = 1 _ ^ > 

- 2(7A)(5A) 

This contradicts our assumption that £1 and £2 are nearly orthogonal. 
So we conclude that pi is close to c, as desired: 



|pi - cl < 8A = 8\/3(n - l)^/5R^/lQl + (Qi/n). (238) 

Finally, using our assumed upper bound on 5, we get that: 

\pi-^< /X. (239) 

So the algorithm succeeds with constant probability. This proves Theorem [6l □ 

Finally, we argue that, when the grid G = (crZ)" H [— L, L)^ is chosen properly, the discrete 
algorithm will behave like the continuous one. Recall from Section 6 that we can approximate a 
function / on R" with a function /a on G, provided that most of the probability mass of / lies within 
distance L of the origin, and most of the probability mass of / lies within distance ^ of the origin. 
This holds for our algorithm, provided that: 

(.„) 

(The first condition follows immediately, since / is supported on a ball. The second condition follows 
from the decay of /; details omitted.) 

Also, we argue that the discretization of the "direction" variable 9 will not introduce too much 
error in the output of the algorithm. The key part of our algorithm involves constructing a line 
£ = {b + \9 \ \ & M}; and if 6 were a continuous variable, this line would pass within distance 
0{^/af3) of the center c. Recall from Section 6 that the discrete curvelet transform resolves 6 within 
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error ib^/a in angular distance. Since b lies at distance 0(/?) from the center c, the error in 9 can 
increase the distance from £ to c by at most 0{^/a|3). Thus the error in the output of the algorithm 
increases by at most a constant factor. 

We can bound the running time of our algorithm as follows. Say we choose L and a so that 
the above inequalities are tight (up to constant factors). Also, suppose 6 satisfies equation (??) 
exactly up to constant factors. Let M be the number of grid points along one direction. Then 
M = 2L/a < 0{B?rc' / ^'^), and the running time is < poly(n,logM) < poly(n, log i?, log ^). 

F.3 Classical lower bound 

We will prove a lower bound on classical randomized algorithms for finding the center of a radial 
function. Recall that an instance of the problem is specified by an oracle / and parameters n, -R, 5 
and fi. Let us define m = lg{R/fi). Consider algorithms that query points within some finite subset 
D C M". We will show that J7(nm/ Ig(nm)) queries are needed to solve this problem. 

Theorem 8 The following holds for any parameters n, R, 5 and /i, and for any finite subset D C M". 
Let us define m = \g{R/ (j). For any classical randomized algorithm that queries points within the 
set D and uses at most (^nm)/ Ig(^nm) oracle queries, there exists a problem instance (depending 
on D, and having the specified parameters n, R, 6 and ^) that causes the algorithm to fail with 
probability at least 1 - 2""™-/^ . 

The assumption that the algorithm queries points within the set D can be understood intuitively 
as follows. We are assuming that the algorithm follows some (arbitrary) convention for how it 
describes points in R" when it queries the oracle. We let D C be the set of points that can be 
queried, and note that D must be finite: if the algorithm runs in time T, then clearly \D\ < 2^. 

Note that this assumption does not weaken our lower bound. The assumption is needed because 
the construction of the hard instance depends on D; however, the actual lower bound (i.e., the 
number of oracle queries and the resulting probability of success) is independent of D. 

Proof: We will use the following version of Yao's minimax lemma, due to Rademacher and Vempala 

m- 

Lemma 9 Let Z be a set of problem instances and A be a set of deterministic algorithms. For any 
probability measure vr over I, and any probability measure v over A, we have that 

inf Pr [algorithm A fails on instance L] < sup Pr [algorithm A fails on instance /]. (241) 

We will use the following approach. First we will fix a probability distribution vr over instances, and 
prove that the best deterministic algorithm still fails with high probability. Note that a randomized 
algorithm is simply a probability distribution over deterministic algorithms; so the minimax lemma 
implies that for any randomized algorithm v, there exists an instance that causes the algorithm to 
fail with high probability. 

We let A be the set of deterministic algorithms that query points within the subset D and make 
at most £ queries, for some i to be specified later. 

We let I be the set of problem instances, such that n, R, 5 and /i are fixed, but / and ccan vary. 
That is, we fix the dimension n, the radius R of the ball in which the center lies, the thickness 5 of 
the spherical shells, and the desired accuracy fi. We define m = lg{R/fi), which is fixed. However, 
the values returned by the radial function /, and the location of its center c, are arbitrary. 
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We now fix a distribution vr on problem instances. Random instances according to this distribu- 
tion are constructed as follows. First, we choose the center point c uniformly at random from the 
ball of radius R around the origin. Let / be of the form 

/(f) = [sf, if X G Ak], where Ak = {x \ kS < \x - c\ < {k + 1)6}, for A; = 0, 1, 2, ... . (242) 

Next, we choose the values of / on the points in the set D (these are the points that the algorithm 
can query). Equivalently, we will choose the values of s^, for those k such that AkD D ^ Let K 
be the set of those k, and note that \K\ < \D\. Let 5 = {1, 2, . . . , |-D|}. Choose a random injective 
map a : K ^ S, and then set Sk = (T{k). 

Consider any deterministic algorithm A G A, and let U be the set of possible outputs of the 
algorithm. We claim that \U\ < 2^'^^. To see this, note that the algorithm A can be described as a 
decision tree, where each node represents a query to the oracle /, and the algorithm chooses which 
branch to follow depending on the oracle's answer. We claim that after seeing the answer to its 
fe'th query, the algorithm can have at most k distinct branches. This is because, while the oracle 
can return many different values, they are meaningless except in cases where the oracle returns the 
same value as it did for a previous query. (Note that, for any permutation a on the range of /, 
we can replace the oracle / with o" o /, to get a new instance of the problem that has the same 
desired solution. These two instances occur with equal probability under the distribution tt.) So, 
when the algorithm receives the answer to its k'th. query, all it can do is compare that value to the 
answers to its previous queries. The number of branches is at most the number of distinct values 
that have been seen previously (which is at most A; — 1), plus 1 (if the new value does not match 
any of the previous ones); thus the number of branches is at most k. Finally, note that the size of 
U is at most the number of leaves at the final level of the decision tree, i.e., after the ^'th query. So 



Thus we can upper-bound the probability that algorithm A succeeds on a random instance 



\U\ < l-2-3---£</ = 2^is^. 




Pr [algorithm A succeeds on instance /] 




volume of ball of radius 



(243) 



Now suppose that £, the number of queries, is at most {^nm) / lg{^nm) . Then 



(Ig(inm) -lglg(i 



nm)) < i 



^nm (assuming nm > 4). 



(244) 




So we have 



Pr [algorithm A succeeds on instance I] < 2 



(245) 



and this holds for all algorithms A G A. So 




(246) 



Now plug into the minimax lemma and the result follows. □ 
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F.4 Finding the center through multiple iterations 

Here we analyze the "continuous" analogue of our algorithm for finding the center of a radial function. 

First, we analyze the procedure OneRound(). Consider what happens for each i £ {1,2}. 

We quantum-sample over a ball of radius R' around 0, then measure the value of /, and get a 
superposition over a shell of radius around c. If the shell is too large, it will not lie completely 
within our original ball, so we only get a fragment of the shell. However, if < R' — R, then we 
are guaranteed to get a complete shell. 

We claim that we observe < Rl - R with probability >l-0{^). To see this, write: 

r ^1 ^1 volume of ball of radius R' — R around c 

Y>t[[3^^> < R - R] = — 

volume of ball of radius R' around (247) 



{R' - RY _ /I J_vn > -(1.2)/S > 1 _ M 



using the fact that 1 — x > e (i-^)^ for all < x < 1/4 (recall that we assumed n > 4). 
Our algorithm does not know the shell's true radius /3^*^ , so it uses 

^(i) =R> -R as an estimate. 
We claim that < < with probability > 1 - O(^). Observe that 

Pr[/3W < {1/S){R! - R)] = '^Q^'^™^ °f ^^^^ ^^^^^^ {^/S){R' - R) around c 

volume of ball of radius R' around (248) 

So 



Pr[(l/5)(i?' -R)< /?» <R!-R] = {1- ( 1 )«)(! - A_) 

>(l-i)(l-M)>l-f. ^'''^ 

Next, we claim that we observe a^*) < ?7*^*\ with constant probability. This follows from Theorem 

II 

Pr[aW < T^W] > 0.045. (250) 
Prom this point on, we take probabilities conditioned on a^*-* < t?*-*^. 

— ~* (i) 

Without loss of generality, assume c = 0. Let H^ be the projector onto the subspace orthogonal 
to 6^^\ and let Hg^ be the projector onto the direction 6^^\ 

We claim that |n^*''5(*)| is smah, and iHg"*?^*-*! is of order /?, with probability > 1 — 0{^). We 
use a similar argument as in the previous section, together with TheoremHl We define e*-*^ = 5/P^^\ 
We get that, with probability > 1 — -I, 



nf)&«| < VsJin - l)e»/?« J(507 + %) • 5, |n^'^6«| < ^23 + %, (251) 



for some constants Qi and Q2- 

Using the definition of e^*^ and the fact that < R' — R = (nS — 1)R < nSR, this imphes 
that 



|nP?»| < S'^/^nV6R^507+Q^, |n^*^6«| < S^I'^nR^2'i + (252) 

The algorithm carries out this procedure twice, for i = 1 and 2. With probability > 1 — O(^), 
this produces two hues, h = {i^^) + A^'^^) | A G M} and £2 = {^^^^ + A^^^) | A G R}, which both pass 
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near the point c. The algorithm then checks that these hnes are nearly orthogonal, and if they are, 
it returns the point on ii closest to ^2- A straightforward calculation shows that this point is given 
by ^^e'W 

First, we claim that the lines £i and £2 are nearly orthogonal {\9^^^ ■ 6^'^^\ < 3/4) with probability 
> 5/9. This follows from the same argument as in the previous section. 

Next, we claim that when ii and £2 are nearly orthogonal (l^^^-* • ^'-^•'l < 3/4), the point on li 
closest to £2 (call it pi) is close to c. Using the same argument as in the previous section, we conclude 
that: 

\pi-^< ^S^I'^nVSR^hm + (253) 
Finally, using our assumed upper bound on (5, we get that: 

|pi - cl < VR^Ji^. (254) 

In summary, we have shown that the procedure OneRound() has the following two properties: 
(1) OneRound() returns a point q (rather than "no answer") with probability 

> (0.045)2 • I > 0.0011, (255) 

and (2) when OneRound() returns a point q, that point q\\es within distance \/S\//u/2 of the center 
point c, with probability 

> ((1-¥)(1- I))' >(1 -¥)'>!-¥• (256) 
We now analyze the complete algorithm, consisting of multiple iterations. 

First, let us upper-bound the number of iterations, assuming that every iteration is successful. 
Let k denote the number of iterations. We claim that k < \\g Ig ^] . 

Let Ri denote the value of Rcur following the z'th iteration, so we have Rq = R (when the 
algorithm starts), Ri+i = ^/Ri\/ /x/2 (the recurrence relation), and Rk < fJ' (when the algorithm 
finishes). 

Let us define Ri = Then we have Rq = (when the algorithm starts), Ri+\ = \fRi (the 
recurrence relation), and R^ < 2 (when the algorithm finishes). 

It is easy to see that R^ = {Rq)^^!'^^^ ^ and a straightforward calculation shows that it suffices to 
set fe= [Iglg^]. 

Next, we show that the algorithm succeeds (i.e., every iteration is successful) with constant prob- 
ability. First, consider a single iteration. The algorithm makes ntries attempts to run OneRound(), 
and it succeeds if at least one of those attempts returns a point q that lies near the center. The prob- 
ability that OneRoundO returns "no answer" every time is < (0.9989)"*"'=» = (0.9989)^^° '"s-^ < 1/5. 
So the probability that OneRound() returns a point q at least once is > 1 — ;5- When this happens, 
the point q is near the center with probability > 1 — ¥• iteration succeeds with overall 

probability > (1 - i)(l - > ^ " ¥ = 1 " 

The algorithm makes n^ter iterations, and it succeeds if all of the iterations succeed. This occurs 
with probability > (1 — 7^;^)"'*'"' > l/e^- (We assumed that R > 8/i, so nuer > 2, and we used the 
fact that 1 — X > e~'^^ for all < 2; < 1/2.) So the algorithm succeeds with constant probability. 

Finally, note that the number of oracle queries is 2ntriesniter < 0(lglg ^ Iglglg ^). This proves 
Theorem [71 □ 
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Finally, we argue that, when the grid G = (crZ)" n [— L, L)" is chosen properly, the discrete 
algorithm will behave like the continuous one. Recall from Section 6 that we can approximate a 
function / on M" with a function /2 on G, provided that most of the probability mass of / lies within 
distance L of the origin, and most of the probability mass of / lies within distance of the origin. 
This holds for our algorithm, provided that: 

H,_L^ (-) 

(The first condition follows immediately, since / is supported on a ball. The second condition follows 
from the decay of /; details omitted.) 

Also, we argue that the discretization of the "direction" variable 9 will not introduce too much 
error in the output of the algorithm. The key part of our algorithm involves constructing a line 
£ = {6 + A0|AgM}; and if were a continuous variable, this line would pass within distance 
0{y/a(3) of the center c. Recall from Section 6 that the discrete curvelet transform resolves within 
error ziz^/a in angular distance. Since b lies at distance 0{f3) from the center c, the error in 9 can 
increase the distance from £ to c by at most 0{-\/al3). Thus the error in the output of the algorithm 
increases by at most a constant factor. 

We can bound the running time of our algorithm as follows. Say we choose L and a so that 
the above inequalities are tight (up to constant factors). Also, suppose 5 satisfies equation (??) 
exactly up to constant factors. Let M be the number of grid points along one direction. Then 
M = 2L/a < 0{B?v? / jjp'), and the running time is < poly(n,logM) < poly(n, log i?, log ^). 
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